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1.0  INTRODUCTION 

Deceleration  measurements  on  a  45-deg  sharp  cone  cylinder  traveling  in  a  constraining 
track  mechanism  at  McDonnell  Douglas  Astronautics  Company  (MDAC)  at  Huntington 
Beach,  California  (Ref.  1),  revealed  increases  in  "apparent"  vehicle  drag.  The  two  obvious 
sources  of  added  drag  appear  to  be  sliding  friction  drag  caused  by  a  portion  of  the  cylinder 
bearing  on  one  or  more  of  the  constraining  tracks  and  drag  caused  by  bow  shock  reflections 
off  the  track  impinging  back  onto  the  model  cone  surface.  The  design  of  future  experiments 
and  tracks  for  the  Hyperballistic  Ranges  (G  and  K)  of  the  von  Karm&n  Gas  Dynamics 
Facility  (VKF)  or  a  full-scale  reentry  facility  depends,  in  part,  on  the  isolation  of  the 
cause  of  this  added  drag  and  the  evaluation  of  the  relative  magnitudes  of  the  aerodynamic 
bow  shock  and  sliding  friction  sources. 

A  simulated  track  mechanism  was  constructed  and  tested  in  Tunnel  F  using  a  sharp 
45-deg  and  a  slightly  blunted  10-deg  cone-cylinder  as  model  configurations.  One-component 
axial-force  measurements  were  made  with  and  without  the  track  mechanism  in  place.  Design 
of  the  track  simulated  range  environment  as  closely  as  possible  and  ensured  that  any 
increase  in  measured  drag  could  be  only  from  sources  induced  by  the  model  bow  shock 
system  and/or  the  shock  system  originating  from  the  rail  leading  edge. 

2.0  TEST  ARTICLE  AND  WIND  TUNNEL 
2.1  TEST  ARTICLE 

Two  cone-cylinder  models  were  tested  with  and  without  a  constraining  track  designed 
from  existing  model  and  track  dimensions  of  the  McDonnell  Douglas  experiment.  Model 
and  track  dimensions  are  given  in  Fig.  1.  A  photograph  of  the  10-deg  model  mounted 
on  the  one-component  force  balance  with  the  track  mechanism  in  place  is  shown  in  Fig. 
2.  Within  constraints,  the  track  rails  could  be  positioned  at  varying  distances  relative  to 
the  force  models.  An  electric  circuit  ensured  that  no  contact  between  model  and  track 
was  made  during  a  test  run.  In  addition  to  these  models,  a  l-in.-diam  hemisphere-cylinder 
with  an  overall  length  of  2.75  in.  was  tested.  This  model  was  instrumented  with  a  stagnation 
point  pressure  gage  which  allowed  local  tunnel  flow  properties  to  be  obtained.  Axial-force 
and  p„  measurements  were  obtained  with  this  model. 

It  was  recognized  that  complete  simulation  of  the  range  environment  was  not  possible 
because  of  the  leading-edge  shock  generated  from  the  track  mechanism  interacting  with 
the  model  and  model  shock.  The  leading  edge  of  the  track  mechanism  was  designed  with 
a  knife  edge  (Figs.  1  and  2)  for  this  reason.  The  influence  of  track  leading-edge  position 
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relative  to  the  model  affected  the  amounts  of  shock  spreading*  and  boundary-layer  growth 
on  the  track  ahead  of  the  model.  These  factors  were  expected  to  influence  the  Tunnel 
F  data  but  would  not  exist  on  the  range  track.  Later  analysis  of  the  data  considers  these 
factors,  and  the  data  are  corrected  for  their  influence. 

2.2  -TUNNEL  AND  NOZZLE  DESCRIPTION 

The  Hypervelocity  Wind  Tunnel  (F)  is  an  arc-driven  wind  tunnel  of  the  hotshot  type 
(Ref.  2)  and  is  capable  of  providing  Mach  numbers  from  about  7.5  to  20  over 
a  Reynolds  number-per-foot  range  from  0.05  x  106  to  70  x  io6.  jest  sections  of  108-in. 
diameter  (M^  =  14  to  20)  and  54-in.  diameter  (M^  =  10  to  17)  are  available  using  a 
4-deg,  half-angle  conical  nozzle.  The  range  of  Mach  numbers  at  a  particular  test  station 
in  the  conical  nozzle  is  obtained  by  using  various  throat  diameters.  The  M^  =  8  and  12 
contoured  nozzles  have  25-in.  and  40-in.  exit  diameters,  respectively,  which  connect  to 
the  54-in.-diam  test  station  and  provide  a  free-jet  exhaust.  The  test  gas  for  aerodynamic 
and  aerothermodynamic  testing  is  nitrogen.  Air  is  used  for  combustion  tests.  The  test 
gas  is  confined  in  either  a  1 .0-ft3 ,  a  2.5-ft3 ,  or  a  4.0-ft3  arc  chamber,  where  it  is  heated 
and  compressed  by  an  electric  arc  discharge.  The  increase  in  pressure  results  in  a  diaphragm 
rupture,  with  the  subsequent  flow  expansion  through  the  nozzle.  Test  times  are  typically 
from  50  to  200  msec.  Shadowgraph  and  schlieren  coverage  are  available  at  both  test 
sections. 

This  test  was  conducted  in  the  1 08-in. -diam  test  section  of  the  conical  nozzle  for 
M,  =  15.4  to  19.2.  Nitrogen  was  the  test  gas.  The  1.0-  and  2.5-ft3  arc  chambers  were 
used,  and  useful  test  times  up  to  approximately  120  msec  were  obtained.  Because  of 
the  relatively  short  test  times,  the  model  wall  temperature  remained  essentially  invariant 
from  the  initial  value  of  approximately  300°K;  thus  Tw/T0  0.1  and  approximates  the 
condition  of  practical  interest  for  reentry  vehicles. 

3.0  INSTRUMENTATION  AND  TEST  PROCEDURES 
3.1  INSTRUMENTATION 

Axial-force  measurements  were  obtained  using  a  one-component,  low  load  internal 
balance  which  is  1 .5  in.  long  and  0.50  in.  in  diameter.  A  calibration  was  performed  before 
the  test  entry  to  determine  whether  interactions  were  present.  None  could  be  detected 
under  the  largest  conceivable  normal  or  side  forces.  Static  loads  from  0.007  to  2.123 


•The  term  "shock  spreading"  refers  to  the  increasing  area  on  the  model  cone  surface  which  is 
influenced  by  the  rail  lip  shocks  as  the  rail  is  extended  forward  of  the  model  shoulder. 


9 


AEDC-TH-75-78 


lbf  were  applied  to  determine  the  balance  calibration  constant.  Absolute  measurement 
uncertainty  was  determined  from  calibration  repeatability  and  is  estimated  as  ±0.001  lbf. 

Test  section  flow  conditions  were  monitored  with  1 .0-in.-diam  hemisphere-cylinders 
instrumented  with  slug  calorimeters  and  pitot  probes.  In  addition  to  the  monitor  probes, 
the  hemisphere-cylinder  model  mentioned  above  was  instrumentated  with  a  sealed, 
miniature  semiconducter  strain  gage-type  pressure  transducer. 

Detailed  information  concerning  the  force,  heat-transfer,  and  pressure  instrumentation 
can  be  found  in  Refs.  3  and  4. 

3.2  PROCEDURES 

The  large  size  of  the  Tunnel  F  test  section  allowed  the  test  to  be  conducted 
concurrently  with  a  sponsored  force  test  using  spare  instrumentation  channels.  The  model 
was  placed  either  13  in.  below  or  12-7/8  in.  above  the  tunnel  centerline,  depending  on 
requirements  dictated  by  the  primary  test  article  located  on  the  nozzle  centerline  (Fig. 
3).  The  model  was  adjusted  to  the  local  velocity  vector,  and  the  model  nose  was  near 
the  window  vertical  centerline  for  all  runs.  Test  variables  included  flow  condition,  model, 
and  track  location.  Table  1  summarizes  the  range  and  combination  of  variables  studied. 


Optional  Viewirc  Ports  -  2  Ports 
at  23  cos  and  1  Part  at  43  deg  Irom 
the  Vertical  ITypl,  Windows  Are 


11  in  in  Diameter 


Optional  Viewing  Ports  -  2  Ports  at  20  deg  and 
2  Ports  at  43  deg  Irom  the  Vertical;  Windows 
Are  11  in  in  Diameter 


4-deg  Conical  Nozzle 
Tnroat  Size  variable 


Optical  Window  Diameter  36  in  ; 
Scltlleren  Picture  30  in.  Diameter 

Access  Door 
Manmum  Area 
Alternate  Optical  Window 
15  5  x  27. 4  in  IMay  Be 
Rotated  in  30-deg  Increments! 


Remote  Angle-ol -Attack 
Adjustment  Irom  -5  to  +20 
deg.  Sector  Angle  ol  Attack 
from  -50  to  +65  deg  by 
Knuckle  Adjustments.  The 
Use  of  Knuckles  Results  In 
Movement  of  the  Model 
Center  of  Rotation  < 


WWW 
Isolation  Mounting 


Figure  3.  Schematic  of  Tunnel  F  108-in.  test  section. 
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Table  1.  Test  Matrix 


Model 

L. 

in. 

(Approximate) 

Re/ft  x  10“6 

4788 

Hemisphere  Cylinder 

Track  Not  Used 

18.5 

0.47  -  0.85 

4789 

Calibration  Body 

19.5 

0.42  -  1.08 

4790 

Calibration  Body 

19.0 

0.39  -  0.65 

4791 

10-deg  Blunt  Cone 

Not  in 

Place 

18.5 

0.49  -  0.80 

4792 

Cylinder 

Not  in 

Place 

19.2 

0.42  -  0.50 

4793 

0. 

50 

19.0 

0.51  -  0.89 

4795 

0. 

75 

19.5 

'0.39  -  0.68 

4812 

0. 

75 

16.2 

0.60  -  1.58 

4816 

0. 

75 

15.5 

0.73  -  1.24 

4800 

1. 

50 

19.0 

0.44  -  0.73 

4794 

3. 

96 

19.0 

0.41  -  0.57 

4799 

3. 

96 

19.2 

0.43  -  0.68 

4796 

3. 

96* 

18.8 

0.39  -  0.70 

45-deg  Sharp  Cone 

Not  in 

Place 

15.4 

0.22  -  0.49 

Cylinder 

15.6 

0.30  -  0.49 

15.6 

0.34  -  0.72 

4805 

1. 

15 

16.0 

0.37  -  0.52 

4804 

0. 

35 

15.6 

0.34  -  0.67 

4803 

0. 

81 

15.8 

0.32  -  0.79 

4808 

2. 

00 

15.3 

0.25  -  0.74 

4819 

2. 

00 

15.7 

0.32  -  0.58 

4807 

3. 

68 

15.5 

0.29  -  0.66 

*Three  Rails  Removed 


3.2.1  Test  Conditions 

The  method  of  determining  the  tunnel  flow  conditions  is  briefly  summarized  as 
follows:  instantaneous  values  of  reservoir  pressure,  p0,  and  free-stream  pitot  pressure,  pj,, 
are  measured,  and  an  instantaneous  value  of  the  stagnation  heat-transfer  rate,  q0,  is  inferred 
from  a  direct  measurement  of  a  shoulder  heat  rate  on  a  l.O-in.-diam  hemisphere-cylinder 
heat  probe.  Total  enthalpy,  Ho,  is  calculated  from  p£,,  q0,  and  the  heat  probe  radius, 
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,  using  Fay-Riddell  theory  (Ref.  5).  The  value  of  Ho  determined  in  this  manner  and  the 
measured  value  of  reservoir  pressure  is  then  used  to  determine  corresponding  values  of 
reservoir  temperature,  density,  and  entropy  from  tabulated  thermodynamic  data  for 
nitrogen  (Ref.  6).  The  reservoir  conditions,  the  measured  value  of  p„,  and  the  assumption 
of  isentropic  flow  in  the  nozzle  are  then  used  to  compute  the  free-stream  conditions. 
The  basic  procedure  followed  in  this  computation  is  given  in  Refs.  7  and  8. 

3.2.2  Data  Acquisition 

Acquisition  of  test  data  was  accomplished  using  a  digital  system  with  concomitant 
oscillographs  for  analog  records.  The  digital  data  system  is  capable  of  scanning  70  data 
channels  in  one  millisecond  and  can  store  up  to  150  scans  of  data.  Basic  data  reduction 
is  done  on  an  off-line  digital  computer. 

4.0  DATA  PRECISION 

The  uncertainties  in  the  calculated  force  data  were  estimated  by  using  the  Taylor 
series  method  of  error  propagation  to  combine  the  uncertainties  in  each  measurement 
occurring  in  the  calculation.  In  general,  it  is  estimated  that  for  nominal  loads  the 
uncertainty  in  the  force  measurement  is  ±6  percent.  This  uncertainty  includes  calibration 
linearity  and  repeatability,  instrumentation  system  error,  and  errors  introduced  by  dynamic 
effects  resulting  from  the  impulsive  operating  nature  of  the  facility.  For  nominal  load 
situations,  the  uncertainty  of  ±6  percent  of  the  balance  component  measurement  combined 
with  a  ±  4-percent  uncertainty  in  the  dynamic  pressure  gives  an  uncertainty  in  Ca(  of 
±7  percent. 

Laboratory  calibrations  using  static  loads  indicate  that  the  pressure  transducers  are 
accurate  to  ±1  percent.  Similarly,  the  uncertainties  in  the  heat-transfer-rate  gages  are  ±5 
percent.  The  uncertainties  in  measured  data,  however,  are  higher  because  of  the  dynamics 
of  the  measurements  and  systems  errors.  The  uncertainties  in  the  monitor  probe 
measurements,  po  and  q0,  and  arc-chamber  measurement,  p0,  were  estimated  considering 
both  the  static  calibrations  and  the  repeatability  of  the  test  section  calibration  profile. 
The  uncertainties  in  the  pressure  data  (po  and  p0)  are  estimated  to  be  ±4  percent  and 
±5  percent,  respectively,  based  on  an  average  of  two  measurements  each,  and  the 
heat-transfer  rate,  q0,  is  ±5  percent  based  on  an  average  of  four  measurements.  These 
values  were  used  to  estimate  uncertainties  in  the  tunnel  flow  parameters  using  the  Taylor 
method  of  error  propagation.  Representative  parameters  are  given  below. 

Uncertainty  (±),  percent 
Re/ft  T,  _p 

1.5  10  6  6 
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The  model  attitude  position  was  set  prior  to  each  run,  and  the  pitch  angles  are 
estimated  to  be  accurate  within  ±0.20  deg. 

The  effect  of  placing  a  model  in  conical  source  flow  was  examined  since  the 
measurements  presented  in  this  report  were  obtained  in  a  conical  nozzle  at  an  off-centerline 
position.  A  source  flow  calculation  indicated  that  a  maximum  error  of  0.6  percent  in 
Ca{  would  result  for  the  10-deg  cone-cylinder  model,  with  smaller  errors  for  the  45-deg 
cone-cylinder  model.  Considering  the  small  size  of  this  correction  and  the  uncertainty 
of  the  force  coefficients  themselves,  these  corrections  were  not  applied  to  the  data 
presented  herein.  All  axial-force  data  represent  total  measured  axial  force.  No  adjustments 
for  base  pressure  effects  have  been  made,  and  no  base  pressure  measurements  were  obtained. 

5.0  RESULTS  AND  DISCUSSION 

Although  it  was  recognized  that  model  bow  shock  location  in  relation  to  the 
constraining  track  wou.ld  be  of  primary  importance,  it  was  not  possible  to  utilize  the 
tunnel  optics  as  an  accurate  measurement  device.  The  demands  of  coverage  on  the  primary 
test  model  located  on  the  tunnel  centerline  and  the  need  for  schlieren  movies  of  this 
model  resulted  in  poor  coverage,  with  few  good  quality  schlieren  stills  of  the  secondary 
models  and  track  system.  No  high  quality  coverage  was  obtained  on  the  45-deg  sharp 
cone-cylinder  model.  Figure  4  indicates  two  of  the  better  shock  pattern  photographs  with 
the  10-deg  cone-cylinder  model  at  two  track  positions.  It  can  be  observed  that  the  model 
bow  shock  was  not  completely  captured  for  the  rail  position  of  0.75  in.  from  the  shoulder, 
whereas  the  1.50-in.  position  exhibited  complete  capture  well  downstream  of  the  rail 
leading  edge.  Although  the  rails  appear  to  have  a  blunt  leading  edge  in  Fig.  4,  this  is 
not  the  case  (see  Figs.  1  and  2).  A  knife  edge  was  fabricated  to  lessen  rail  leading-edge 
shock  effects  since  a  true  range  environment  would  not  produce  such  a  shock  system. 
It  will  be  demonstrated  in  Section  6.0  that  the  schlieren  coverage  shown  in  Fig.  4  was 
very  important  in  analysis  of  the  data  and  its  application  to  a  full-scale  track  system. 

Axial-force  measurements  on  the  hemisphere-cylinder  calibration  body  were  in 
excellent  agreement  with  previously  published  data  (Ref.  9).  These  data  are  not  presented 
herein  since  they  would  serve  no  useful  purpose  in  the  present  analysis. 

In  order  to  define  a  "baseline"  value  for  later  determination  of  the  effect  of  the 
track  rail  system,  several  runs  were  made  with  the  10-deg  blunt  cone-cylinder  and  the 
45-deg  sharp  cone-cylinder  with  the  track  mechanism  completely  removed.  Typical  data 
from  these  runs  are  tabulated  in  Table  2  and  shown  in  Fig.  5.  The  dimension  used 
in  calculation  of  the  parameter  M J\f  Reg  is  defined  for  each  model  in  this  figure. 
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Run  4800,  L  =  1.50  in 


*  — — —  «* 

Figure  4.  Schlieren  coverage  of  the  10-deg  cone-cylinder  model 
with  two  rail  positions. 
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Table  2.  Flow  Conditions  and  Axial-Force  Data  with 
No  Track  System  in  Place 


■ 

Time, 

■gn 

T 

It/ sec 

Re/ ft 

c*t 

Model 

E 

msec 

psla 

°r’ 

x  10~6 

10-deg  Blunt 

4791 

50 

22.764 

4,667 

0.00525 

75.4 

8,324 

1.357 

0  .799 

0.126 

60 

20,758 

5,121 

0.00500 

85.9 

8,716 

1.246 

18  .9 

0.615 

0.135 

80 

17,061 

4,885 

0.004  52 

83.1 

8,  467 

1.099 

18.7 

0.576 

0.152 

90 

15,624 

4,837 

0.00436 

83.3 

8,  407 

1.043 

18.5 

0.551 

0.159 

ins 

14,109 

4,725 

0.00420 

82.4 

8,289 

0.987 

18.3 

0.534 

0.163 

110 

13,010 

4,759 

0.00396 

83.6 

8,307 

0.921 

18.2 

0.490 

0.166 

4792 

90 

15,443 

4  775 

0.00333 

76.1 

8,353 

0.860 

19.2 

0.500 

0.152 

4792 

110 

12.728 

4.730 

0.00272 

75.0 

8.284 

0.700 

19.2 

0.417 

0.152 

45-de{ 

Cone-C; 

Sharp 

/Under 

48 

02 

60 

70 

5,363 

4,909 

0.00494 

0.00463 

88.0 

103.9 

7,316 

7,837 

0.848 

0.771 

mu 

0.486 

0.350 

1 .066 
1.055 

80 

4,56  5 

4,810 

0.00427 

116.7 

8,245 

0.700 

15.3 

0.269 

1.083 

90 

4,177 

4,858 

0.00382 

117.4 

8,286 

0.628 

15.3 

0.239 

1.092 

100 

3,860 

4,689 

0.00338 

110.9 

8,127 

0.566 

15  5 

0.232 

1.083 

110 

3,552 

4,570 

0.00305 

106.7 

8,011 

0.515 

1  5.6 

0.223 

1.075 

4806 

80 

4, 529 

3,548 

0.00433 

6.986 

0.741 

15.6 

0.486 

0.961 

100 

3,880 

3,897 

0.00365 

89.8 

7,344 

0.618 

15.6 

0.346 

0.988 

110 

3,  572 

3,901 

0.00339 

Blfll 

7,344 

0.572 

15.5 

0.319 

0.961 

120 

3,329 

3,901 

0.00300 

88.8 

7,343 

0.513 

15.6 

0.291 

0.973 

4817 

70 

5,439 

3,209 

0.00529 

71.2 

6,622 

0.916 

15.7 

0.718 

0.928 

80 

5,084 

3,638 

0.00478 

82.4 

7,087 

0.820 

ESI 

0.519 

0.934 

90 

4,669 

3,8  57 

0.00431 

88.2 

7,313 

0.737 

0.414 

0.935 

100 

4,340 

4,034 

0.00391 

92.5 

7,490 

0.667 

B31 

0.356 

0.933  | 

110 

4,011 

3,944 

0.00357 

89.7 

7,395 

0.613 

15.7 

0.342 

0.935  | 

Since  the  present  measurements  were  made  in  relatively  low  free-stream  Reynolds 
number  conditions,  second-order  viscous  shear  stress  effects  on  the  aerodynamic  drag 
influenced  the  results.  The  author  could  find  no  previous  data  on  blunt  cone-cylinders 
in  this  particular  flow  regime  although  there  are  numerous  sets  of  data  available  on  blunt 
cones.  Since  the  cylinder  portion  of  the  10-deg  blunt  cone  model  represented  about  25 
percent  of  its  axial  length,  it  could  be  expected  that  this  portion  of  the  model  would 
increase  the  viscous  drag  significantly.  Previous  data  on  blunt  cones  suggest  that  total 
drag  should  be  from  1.4  to  1.7  times  the  inviscid  pressure  drag.  A  calculation  of  this 
latter  parameter  (Ca*)  using  the  method  of  Ref.  10  and  a  Newtonian  proportionality 
constant  of  1 .833  gave  a  value  of  0.081 .  Thus,  based  on  previous  cone  data  and  approximate 
analytic  results,  a  measured  axial-force  coefficient  of  0.11  to  0.14  could  be  expected  on 
the  10-deg  blunt  cone  cylinder  with  no  track  interference.  The  influence  of  viscous  drag 
on  the  cylindrical  portion  of  the  model  would  increase  this  value  somewhat. 

For  the  case  of  the  45-deg  sharp  cone-cylinder  model,  previous  data  from  the  VKF 
Low  Density  Hypersonic  Wind  Tunnel  (L)  and  Tunnel  C  and  a  Langley  hypersonic  tunnel 
using  40-,  45-,  and  50-deg  sharp  cones  suggest  that  very  little  if  any  viscous  contribution 
to  the  total  drag  would  be  present  at  the  free-stream  Reynolds  numbers  of  the  present 
test.  An  inviscid  solution  from  Ref.  11  indicates  that  an  axial-force  coefficient,  Cap  of 
1.06  would  be  measured  on  a  45-deg  sharp  cone  with  no  cylindrical  afterbody.  The  viscous 
drag  on  the  cylindrical  afterbody  could  be  expected  to  increase  this  value  slightly. 
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Data  from  Runs  4802,  4806,  and  4817  (Table  21 
4*  4j  +  l 2  -  2.4°0  in. 


X 

-  h 

45-deg  Sharp  Cone-Cylinder 


Data  from  Runs  4791  and  4792  (Table  2} 


The  data  for  both  the  10-deg  slightly  blunted  cone-cylinder  and  the  45-deg  sharp 
cone-cylinder  are  in  good  agreement  with  these  analytic  and  experimental  estimates.  The 
fairings  of  the  experimental  results  shown  in  Fig.  5  serve  as  the  baseline  values  for  later 
analysis  of  the  effect  of  the  simulated  track  mechanism  on  shock-induced  (shock 
interference)  drag. 

Results  of  axial-force  measurements  with  the  45-deg  sharp  cone-cylinder  model  with 
the  track  system  in  place  at  varying  locations  relative  to  the  model  shoulder  are  tabulated 
in  Table  3  and  shown  in  Fig.  6.  Although  viscous  drag  influence  was  small,  the  data 
are  shown  as  a  function  of  M JV  Reg  for  convenience. 
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Table  3.  Flow  Conditions  and  Axial-Force  Data  on 
45-deg  Sharp  Cone-Cylinder  with  Track 


Run 

Po, 

psla 

To, 

OR 

Poo  f 
psia 

Too  , 
OR 

Voo 

’  J 

ft/sec 

Qoo  , 

psia 

Ko 

Re/ft 
x  10-6 

E8 

in 

4805 

80 

4,812 

3, 

476 

0.00419 

76.1 

6, 

911 

0.739 

15.9 

0.519 

0. 

15 

1.027 

90 

4,426 

3, 

352 

0.00379 

72.5 

6, 

772 

0.675 

16.0 

0.508 

1.019 

100 

4,122 

3, 

249 

0.00354 

69.8 

6, 

653 

0.632 

16.0 

0 . 503 

1.019 

110 

3,826 

3, 

326 

0.00345 

72.9 

.6, 

735 

0.605 

15.8 

0.455 

1.024 

120 

3,559 

3, 

645 

0.00349 

83.6 

7. 

077 

0.588 

15.5 

0.367 

* 

1.028 

4804 

60 

5,571 

3, 

362 

0.00581 

77.0 

6, 

793 

0.981 

15.5 

0.670 

0. 

35 

1.167 

70 

5,155 

3, 

807 

0.00523 

89.1 

7, 

264 

0.872 

15.4 

9.498 

1.161 

80 

4,788 

4, 

126 

0.00464 

97.2 

7, 

587 

0.774 

15.4 

0.388 

1.185 

90 

4,413 

4, 

100 

0.00408 

95.2 

7, 

559 

0.689 

15.  5 

0.354 

1.192 

100 

4,071 

3, 

967 

0.00363 

90.4 

7, 

420 

0.622 

15.7 

0.343 

1.171 

4803 

70 

4,687 

2, 

812 

0.00453 

60.6 

6, 

152 

0.797 

15.9 

0.791 

0. 

81 

1.244 

80 

4,406 

3, 

254 

0.00411 

71.7 

6, 

662 

0.717 

15.8 

0.585 

1.252 

90 

4,103 

3, 

483 

0.00367 

77.0 

6, 

911 

0.642 

15.8 

0.446 

1.254 

100 

3,823 

3, 

649 

0.00329 

80.7 

7, 

087 

0.577 

15.8 

0.373 

1.254 

110 

3,542 

3, 

773 

0.00302 

83.8 

7, 

213 

0.528 

15.8 

0.323 

1 

1.257 

4808 

70 

4,817 

3, 

137 

0.00532 

71.9 

6, 

531 

0.889 

15.5 

0.735 

2.00 

1.139 

80 

4,542 

3, 

753 

0.00479 

88.5 

7, 

200 

0.790 

15.4 

0.459 

1.116 

90 

4,172 

4, 

073 

0.00434 

97  .6 

7, 

525 

0.709 

15.3 

0.357 

1.148 

100 

3,853 

4, 

309 

0.00400 

104.7 

7, 

757 

0.647 

15.2 

0.294 

1.167 

110 

3,572 

4, 

495 

0.00378 

111.0 

7, 

933 

0.604 

15.1 

0.254 

1.176 

4819 

70 

5,442 

3, 

541 

0.00506 

79.4 

6, 

988 

0.876 

15.7 

0.584 

2.00 

1.148 

80 

5,066 

4, 

059 

0.00449 

92.8 

7, 

525 

0.770 

15.7 

0.415 

1.158 

90 

4,647 

4, 

135 

0.00398 

94.1 

7, 

597 

0.688 

15.7 

0.356 

1.157 

100 

4,279 

4, 

018 

0.00358 

90.2 

7, 

474 

0.625 

15.8 

0.343 

1.147 

T 

iio 

3,964 

4, 

043 

0.00334 

90.5 

7, 

475 

0.581 

15.8 

0.318 

f 

1.137 

4807 

60 

5,117 

3, 

334 

0.00569 

77.6 

6, 

755 

0.942 

15.4 

0.664 

3.68 

1.229 

70 

4,761 

3, 

796 

0.00494 

89.4 

7, 

249 

0.818 

15.4 

0.467 

1.227 

80 

4,454 

4, 

018 

0.00429 

93.8 

7, 

473 

0.719 

15.5 

0.388 

1.253 

90 

4,093 

4, 

149 

0.00381 

96.8 

7, 

604 

0.641 

15.  5 

0.322 

1.294 

100 

3,813 

4, 

126 

0.00343 

95.1 

.7, 

577 

0. 583 

15.6 

0.298 

1.297 

110 

3,519 

4, 

027 

0.00318 

92.5 

7, 

474 

0.542 

15.6 

0.290 

1.259 

A  shadowgraph  supplied  by  McDonnell  Douglas  on  the  45-deg  models  indicated  that 
the  interaction  of  the  bow  shock  with  the  track  would  be  a  Mach  reflection  and  that 
a  triple  intersection  would  occur.  A  sketch  of  the  shadowgraph  is  shown  in  Fig.  7.  Scaling 
this  figure  produced  a  distance  of  0.35  in.  from  the  shoulder  for  the  intersection  of  the 
Mach  reflection  and  track.  An  examination  of  the  schlieren  results  from  run  4805  (L 
=  0.15)  indicated  that  no  interaction  of  the  model  bow  shock  and  the  "inside"  of  the 
rail  occurred.  The  small  increased  drag  for  small  value  of  L  is  therefore  not  surprising. 
The  quality  of  the  schlieren  from  run  4804  (L  =  0.35)  did  not  allow  a  good  determination 
of  shock  interaction  location.  However,  Fig.  7  (from  lower  Mach  number  flow)  indicated 
that  the  interaction  should  have  been  completed  at  this  track  position.  A  17-percent 
increase  in  aerodynamic  drag  was  measured  from  this  run.  Track  data  from  McDonnell 
Douglas  on  a  45-deg  sharp  cone-cylinder  indicated  approximately  a  26-percent  increase, 
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but  friction  forces  between  the  track  rails-  and  the  cylinder  portion  of  the  model  may 
have  been  present  in  the  track  experiments,  and  flow  conditions  (M^  and  Re/ft)  were 
considerably  different. 


Ma>/VR®£ 


Figure  6.  Axial-force  measurements  on  the  45-deg  sharp 
cone-cylinder  with  track  system  in  place. 

The  observed  increase  in  total  axial  drag  in  Fig.  6  as  the  track  leading  edge  was 
further  extended  suggests  a  strong  influence  of  the  shock  system  induced  by  the  rail  leading 
edge  and  boundary-layer  growth  on  the  rail  surface.  Neither  of  these  would  be  present 
in  an  actual  aeroballistic  range  environment,  and  they  would  prevent  direct  application 
of  the  present  data  to  a  range  track  without  modification. 

Results  of  axial-force  measurements  with  the  10-deg  blunt  cone-cylinder  model  with 
the  track  system  in  place  at  varying  locations  relative  to  the  model  shoulder  are  tabulated 
in  Table  4  and  shown  in  Fig.  8.  Since  this  model  was  influenced  by  viscous  shear  forces, 
the  data  are  plotted  as  a  function  of  M JyJ Reg. 
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Figure  7.  Sketch  of  shadowgraph  from  McDonnell  Douglas  track  test. 


Table  4.  Flow  Conditions  and  Axial-Force  Data  on 
10-deg  Blunt  Cone-Cylinder  with  Track 


Run 

Time, 

msec 

Po» 

psia 

To, 

»R 

P»t 

psia 

°R 

V„, 

f  t/sec 

psia 

Re/ft 
x  10“b 

L, 

in . 

cAt 

1793 

60 

19,924 

3, 

985 

0.00390 

59 .0 

7, 

640 

1 .087 

20.0 

0.891 

0. 

50 

0.140 

70 

18,542 

4, 

172 

0.00385 

63.6 

7, 

811 

1.041 

19.7 

0.775 

0.144 

80 

16,725 

4, 

358 

0.00395 

69.6 

7, 

970 

1  .014 

19.2 

0.676 

0.146 

100 

15,897 

4, 

603 

0.00404 

79.1 

8, 

172 

0.960 

18.4 

0.549 

0.143 

110 

12,529 

4, 

586 

0.00381 

79.8 

«, 

142 

0.892 

18.3 

0.507 

0.147 

4795 

80 

16,58  5 

4, 

077 

0.00296 

59.1 

7, 

697 

0.836 

20.1 

0.679 

0.75 

0.176 

90 

15,339 

4. 

352 

0.00294 

65.  5 

7, 

955 

0.799 

19.7 

0.567 

0.183 

100 

14,105 

4, 

646 

0.00293 

72.7 

8, 

220 

0.766 

19.3 

0.473 

0,187 

110 

12,845 
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Figure  8.  Axial-force  measurements  on  the  10-deg  blunt 
cone-cylinder  with  track  system  in  place. 


The  increase  in  aerodynamic  drag  due  to  the  reflection  of  the  model  bow  shock 
from  the  rail  back  onto  the  model  is  seen  to  be  quite  significant!  From  an  examination 
of  the  schlieren  results,  it  was  concluded  that  the  bow  shock  was  not  completely  captured 
for  L  =  0.750  in.  and  that  it  interacted  just  downstream  of  the  track  leading  edge  at 
•L  =  1.50  in.  (Fig.4).  The  fact  that  the  drag  continued  to  rise  as  the  leading  edge  of  the 
track  extended  forward  suggested  again  that  the  track  leading-edge  shock  and/or  the 
boundary-layer  growth  on  the  track  at  the  location  of  the  model  did  violate  simulation 
of  the  range  environment.  In  this  regard,  these  data  are  discussed  in  more  detail  in  Section 
6.0.  There  are  as  yet  no  range  data  to  compare  to  these  measurements. 

6.0  APPLICATION  OF  TUNNEL  F  DATA  IN  DETERMINATION  OF  SHOCK 
INTERFERENCE  DRAG  ON  A  FULL-SCALE  TRACK  SYSTEM 

The  total  drag  coefficient,  Cdt,  experienced  by  a  spherically  blunted  cone-cylinder 
launched  at  zero  angle  of  attack  in  a  range  and  constrained  by  a  track  system  is  the 
sum  of  the  component  drag  coefficients  expressed  as 
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aerodynamic  forces 


track  related  forces 


Cdt  C-Dp  Coy  +  Cq  g  +  C[)£p  +  +  CD] 


_  Total  Drag 
q  A 

*00 


(1) 


The  first  four  components  are  standard  aerodynamic  forces-namely  pressure  drag, 
viscous  drag,  base  drag,  and  viscous  induced  pressure  drag.  The  term  Cdsi  is  the 
contribution  caused  by  the  influence  of  the  upstream  shock  system  on  the  track  rail  and 
the  subsequent  interaction  on  the  model  surface.  The  term  Cdf  is  the  contribution  of 
the  sliding  friction  force  between  the  cylinder  and  the  rails  of  the  track.  The  present 
experimental  program  attempted  to  measure  the  term  Cdsi  for  two  particular  cone 
geometries.  The  term  Cdf  was  intentionally  kept  at  zero. 


It  was  proposed  that  the  Tunnel  F  experimental  data  be  scaled  to  other  track 
dimensions  and  other  model  geometries  and  a  parametric  study  be  conducted  to  determine 
the  magnitude  of  the  shock  interaction  drag  at  typical  range  conditions.  It  was  also 
requested  that  the  study  be  limited  to  sharp  or  slightly  blunted  slender  cones. 


Each  of  the  terms  in  Eq.  (1)  will  be  discussed.  However,  the  method  by  which  the 
term  Cdsi  is  determined  will  be  discussed  in  some  detail.  Future  studies,  with  additional 
experimental  data,  will  no  doubt  improve  the  present  analysis.  Terms  related  to 
track-induced  drag  forces  are  discussed  first. 


6.1  SHOCK  INTERFERENCE  DRAG 


The  evaluation  of  the  term  Cdsi  in  Eq.  (1)  and  its  relative  importance  to  the  other 
terms  of  the  equation  is  the  primary  purpose  of  the  present  study.  The  basic  mechanism 
of  the  induced  drag  on  the  model  due  to  the  reflection  of  the  model  bow  shock  from 
the  rail  back  on  to  the  surface  of  the  model  and  subsequent  reflections  between  the  model 
and  the  track  rail  can  be  seen  in  Fig.  9.  This  analysis  and  discussion  is  limited  to  sharp 
or  slightly  blunted  slender  cones  at  Mach  numbers  sufficiently  high  that  the  local  flow 
in  the  vicinity  of  the  track  rail  is  supersonic  and  the  reflection  is  of  a  "like  sense"  (Ref. 
13)  so  that  the  angle  of  reflection  equals  the  angle  of  incidence.  This  stipulation  removes 
the  experimental  data  obtained  using  the  45-deg  sharp  cone  from  the  present  discussion. 
The  ranges  of  variables  suggested  for  study  were  5  <  0e  <  30  deg  and  0  <  £  <  0.3. 
Therefore,  the  analysis  does  depend  greatly  on  the  data  obtained  using  the  slightly  blunted 
(£  =  0.167)  10-deg  cone  model.  For  a  model  of  finite  bluntness,  the  shock  approaching 
the  track  rail  is  neither  straight  nor  two-dimensional.  A  rigorous  analysis  should  account 
for  shock  curvature  and  local  flow  gradients.  An  example  is  shown  in  Fig.  10  which 
compares  perfect  gas  method  of  characteristics  solutions**  (M^  =  14  and  20)  to  inviscid 


♦♦Solutions  obtained  by  E.  R.  March'and,  Special  Studies  Group,  VKF. 
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flow-field  sharp-cone  solutions  (Ref.  11).  This  figure  is  drawn  to  scale  for  the  present 
10-deg  blunt  cone  with  the  track  leading  edge  placed  1.50  in.  from  the  cone  shoulder. 
A  previously  discussed  schlierien  photograph  (Fig.  4)  obtained  during  run  4800  indicates 
good  agreement  between  observed  and  calculated  shock  shape. 

Side  and  Bottom  Rails  Are  Not  Shown 


Figure  10.  Comparison  of  blunt  and  sharp  cone  shock  shapes 
in  reflection  region. 

Although  it  was  recognized  that  a  blunt  bow  shock  analysis  would  be  of  more  practical 
value,  it  was  not  attempted  because  a  large  number  of  characteristic  solutions  would  be 
required.  Also,  shock  shape  is  sensitive  to  real  gas  effects,  but  all  calculations  in  the  present 
analysis  assume  perfect  gas  relationships.  An  example  of  perfect  and  real  gas  shock  shapes 
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is  shown  in  Fig.  11.  It  was  felt  that  the  lack  of  experimental  data  would  prevent  an 
evaluation  of  the  rigorous  approach  even  if  the  increased  effort  was  undertaken  to  obtain 
the  numerous  necessary  solutions.  The  present  analysis  assumes  that  the  local  pressure, 
ps  (Fig.  9),  and  the  geometric  length,  x,  can  be  computed  by  crossing  two  inviscid  flow-field 
conical  shock  systems.  It  should  be  noted  that  the  reflected  shock  would  more  nearly 
be  plane  since  it  is  from  a  two-dimensional  surface.  However,  using  this  model  greatly 
simplifies  the  numerical  work  required  for  the  analysis. 


Figure  11.  Perfect  and  real  gas  shock  shapes. 


As  Fig.  9  indicates,  the  drag  (neglecting  the  lip  shock)  induced  by  the  bow  shock 
reflection  system  can  be  written 

DSi  =  4K  (x)  (ps/pj  (pj  (RW)  (2) 

and 

Cdsi  =  Dsi/q^A  (3) 


where 

(Ps/P.)(Pj  *s  the  average  pressure  in  the  reflection  region  and  is  assumed 
to  be  the  value  after  the  first  reflection  (subsequent  reflections  ignored); 
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2.  x  is  the  cone  surface  distance  from  the  first  reflection  to  the  shoulder  (Fig. 

9),  RW  is  the  width  of  the  rail,  and  the  constant  4  is  included  because 
a  four  rail  track  system  was  used  (The  term  4(x)(RW)  is  therefore  the 
affected  shock  area  influenced  by  the  pressure  p5); 

3.  K  is  a  complex  correction  factor  which  must  be  empirically  derived  from 
available  data  to  account  for  the  simplification  of  the  model,  lip  shock 
contribution,  and  second-order  effects.  (The  parameter  should  also  allow 
a  range  of  model  configurations  and  track  dimension  to  be  studied.) 

This  model  implicitly  assumes  a  step  change  in  local  surface  pressure,  ps,  down  to 
a  normal  cone  surface  pressure,  pw ,  at  the  edge  of  each  rail  (Fig.  9).  It  also  ignores  the 
influence  caused  by  the  model  bow  shock  interacting  with  the  sides  of  the  rail  and  thereby 
causing  a  local  increase  in  model  pressure.  The  empirical  constant  K  could  be  expected 
to  be  a  strong  function  of  the  product  (x)(ps /p^),  cone  angle  and  bluntness,  and  Mach 
number,  and  a  weaker  function  of  rail  width  and  rail  depth.  Experimental  data  presented 
earlier  served  as  the  basis  for  evaluating  the  constant  K  in  Eq.  (2).  Since  this  is  basic 
to  the  present  analysis,  a  discussion  of  its  evaluation  is  given  in  some  detail. 

To  be  completely  valid;  this  constant  should  be  used  only  for  the  geometry  and 
flow  conditions  for  which  the  data  were  obtained;  i.e., 

15.2  <  <  20.1 

0C  =  10  deg 
£  =  0.167 

laminar  boundary-layer  conditions 

The  present  task  demanded  that  a  method  be  devised  which  would  allow  calculation 
of  Dgi  for  other  geometric  and  flow  conditions;  i.e., 

5  <  0C  <  30 
6.5  "<  M  <20 

oo 

0  <  £  <  0.3 

turbulent  boundary-layer  conditions  and 
range  temperatures  and  pressure  conditions 

The  method  which  is  used  in  accomplishing  this  is  also  discussed  below.  The  need 
for  additional  experimental  data  will  be  apparent  in  this  regard. 
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6.1.1  Determination  of  the  Constant  K  * 

The  drag  measurements  discussed  in  Section  5  and  shown  in  Fig.  8  are  used  in  the 
determination  of  the  parameter  K.  Intuitive  reasoning  suggests  that  the  track  location  which 
"more  nearly"  simulates  a  range  environment  would  be  at  the  position  for  which  the 
model  bow  shock  just  intersects  the  bottom  of  the  track  rail.  In  other  words,  complete 
capture  of  the  shock  system  is  desired  with  no  excess  track  length  to  allow  bounday-layer 
growth  or  lip  shock  effects  to  excessively  influence  the  data.  Three  runs  were  made  with 
the  track  leading  edge  0.750  in.  in  front  of  the  cone  shoulder,  and  a  consistent  20-percent 
increase  in  drag  was  measured.  However,  a  post-test  analysis  of  both  characteristic  blunt 
body  shock  shape  solutions  and  schlierien  photographs  strongly  suggests  that  incomplete 
capture  of  the  shock  system  resulted  at  this  rail  location  (Figs.  4  and  10). 

A  cross  plot  of  the  10-deg  cone  drag  measurements  as  a  function  of  rail  position 
is  shown  in  Fig.  12.  From  characteristics  solutions  (Fig.  10)  it  was  inferred  that  placing 


‘  0  2  4 

L,  in. 

Figure  12.  Interpretation  of  10-deg  cone  data  used 
in  solution  of  the  parameter  K. 
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the  rail  at  about  1.15  in.  upstream  of  the  cone-cylinder  function  would  have  been  the 
optimum  location.  The  manner  in  which  the  data  are  presented  in  Fig.  12  implicitly  assumes 
that  the  observed  increase  in  drag  which  resulted  when  the  rails  were  in  position  was 
entirely  due  to  the  increased  pressure  on  the  cone  induced  by  the  model  bow  and  rail 
lip  shock  and  that  none  of  the  other  terms  in  Eq.  (1)  was  affected  by  the  presence  of 
the  track.  Therefore,  the  differential  between  the  baseline  data  and  the  total  drag  would 
represent  the  contribution  of  of  the  track  system  and  be,  by  definition.  CDs].  From  Fig. 
12,  at  L  =  1.15,  this  results  in 

Cdsi  -  1.34  Cqbl  Cdbl  =  0-34  Cdbl  (4) 

where  CDbl  can  be  read,  for  a  given  value  of  M Jy/ Reg,  from  Fig.  5.  However,  the 
contribution  of  the  rail  lip  shock  is  also  included  in  the  total  shock  interaction  drag 
component,  and  it  would  therefore  be  conservative  to  apply  this  result  to  an  actual  track 
system.  To  correct  for  this  tunnel-induced  effect,  it  was  assumed  that  once  the  bow  shock 
was  completely  captured,  no  further  increase  in  drag  due  to  bow  shock  interaction  would 
result,  t  However,  lip  shock  drag  would  continue  to  increase  with  L  since  boundary-layer 
and  shock  spreading  would  increase.  By  making  the  additional  assumption  that  the  increase 
over  the  drag  measured  at  L  =  0.50  was  due  largely  to  rail  lip  shock  drag,  one  could 
estimate  the  individual  contributions  of  model  bow  shock  and  lip  shock.  The  technique 
is  shown  graphically  in  Fig.  1 2.  In  the  range  from  1 .20  to  1 .34,  various  values  of 

^°BL  +  ^Dfiow  SHOCK^DBL 

were  tried  until  a  reasonable  variation  of  lip  shock  drag  with  rail  position  L  was  found. 
The  resulting  value  gave  a  correction  of  0.09  from  the  total  measured  value  of  1.34,  which 
is  about  a  7-percent  correction  at  L  =  1.15  in.  Although  the  technique  and  assumptions 
used  are  not  by  any  means  rigorous,  the  resulting  correction  is  not  large,  and  a  better 
method  could  probably  not  be  devised  with  the  given  experimental  data.  A  final  value 
of  CDsi,  where  only  the  bow  shock  is  considered,  is  therefore 

Cdsi  =  °-25  CdBl  (5) 

and 

Dsi  =  (0.25  CDbl)  (q^A)  (6) 


+  This  interpretation  of  the  data  was  suggested  by  Mr.  Glen  Norfleet,  Manager  of  the  VKF 
Aeroballistic  Branch. 
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For  a  given  Tunnel  F  run,  where  Cdbl  is  known,  the  parameter  K  can  be  obtained 
from  Eqs.  (2),  (3),  and  (6)  and  is  expressed  as 

0.25  CDbl 

K  =  -  (7) 

4(Ps/pJ(P  J  (x)  (RW) 

Equation  (7)  was  evaluated  at  discrete  time  points  for  runs  4795,  4800,  and  4816, 
which  represented  the  extreme  range  of  q^,  M^,  and  p^  used  during  the  present  test. 
The  term  pjp^  and  x  were  hand  calculated  from  the  conical  inviscid  flow-field  solutions 
of  Jones  (Ref.  11).  The  results  were  as  follows: 

Number  of  . 


Run 

Time  Points 

Average  Mo 

Average  K 

4795 

5 

19.42 

0.52  ±  0.03 

4800 

6 

18.90 

0.49  ±  0.04 

4816 

4 

15.45 

0.35  ±  0.03 

This  variation  of  K  was  initially  very  disturbing  since  the  observed  strong  Mach  number 
influence  was  not  expected.  A  value  for  K  of  unity  would  have  resulted  if  all  flow  model 
assumptions  were  correct  and  second-order  effects  such  as  interaction  of  the  shock  off 
the  rail  side  wall,  boundary-layer  change,  base  pressure  effect,  etc.  were  small.  Considering 
the  relative  simplicity  of  the  flow  model,  the  absolute  values  from  0.35  to  0.52  were 
not  surprising. 

The  variation  of  K  with  Mach  number  was  explained  when  it  was  realized  that  K 
was  very  sensitive  to  changes  in  the  product  (xXPs/P.)-  A  two-step  conical  shock  solution 
was  performed  for  a  10-deg  sharp  cone  in  the  range  3  <!  <^30,*  and  the  resulting 

variation  of  (x)(ps/poo)  is  shown  in  Fig.  13,  as  is  the  variation  of  K.  Dimensional  analysis 
suggested  use  of  the  parameter  M^  sin  0C  as  the  dependent  variable.  Since  the  calculation 
was  performed  for  a  specific  cone  geometry,  the  parameter  (xXPs/p,^  has  the  units  of 
length.  However,  only  the  slope  and  not  the  absolute  value  is  of  interest,  and  the  results 
could  therefore  be  used  for  other  cone  base  diameters.  The  parameter  K,  in  part,  corrects 
the  flow  model  for  the  error  induced,  by  using  sharp  conical  shock  relationships  rather 
than  blunt  body  characteristics  solutions  in  calculating  the  shock  area  4(x)(RW)  and  local 
pressure,  ps.  One  would  expect  this  error  to  increase  as  Mach  number  decreases,  as  is 
indicated  by  the  results  in  Fig.  13.  A  few  experimental  points  at  low  Mach  numbers  (Mm 
«  8)  would  verify  this.  The  experimentally  derived  values  of  K  plotted  as  a  function 
of  Mm  sin  6C  were  then  faired  using  the  slope  of  the  (x)(Ps/p„)  curve.  A  curve  fit  of 
the  variation  of  K  shown  in  Fig.  13  yields  the  empirical  cubic  relationship 


tThe  wide  Mach  number  variation  was  necessary  for  reasons  explained  later. 
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K  w  0.07795  -  0.02892  sin  9C  +  0.05204  sin*  9C  -  0.000841  Mj  sin3  9C  (8) 

It  should  be  noted  that  this  relationship  is  valid  only  in  the  range  3  <  M^  <  30 
and  for  a  10-deg  blunt  cone  with  £  =  0.167. 


Figure  13.  Empirical  constant  K  as  determined  from  Tunnel  F  data 
and  conical  inviscid  flow-field  solutions. 

6.1.2  Derivation  of  K  for  Cone  Angles  Other  Than  10  deg  with  £  =  0.167 

Since  it  is  obvious  that  the  term  (x)(ps/poo)  is  sensitive  to  Mach  number  for  a  given 
cone  angle,  it  was  assumed  that  for  a  given  Mach  number  it  would  be  sensitive  to  cone 
angle.  An  evaluation  of  shock  interference  drag  on  models  with  cone  angles  between  5 
and  30  deg  was  desired.  To  permit  this,  two-step  inviscid  flow-field  sharp  cone  solutions 
using  the  tables  of  Ref.  1 1  were  again  performed  (for  6  <  M^  <  20),  and  the  results 
are  shown  in  Fig.  14.  It  can  be  observed  that  the  ratio  of  (x)(ps/poo)  for  a  cone  angle 
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0C  to  that  for  a  10-deg  cone  is  independent  of  sin  0C  (Fig.  15).  If  it  is  then  again 
assumed  that  the  empirical  constant  K  is  a  strong  function  of  (xKps/p^)  and  Mach  number 
and  only  a  weak  function  of  track  width,  Eq.  (8)  can  be  corrected  for  cone  angles  other 
than  10  deg  by  the  empirical  parameter  Y  defined  in  Fig.  15  and  given  by  the  relationship 

Y  •  8. 1775/(0C)°  9105  (9) 


where  0C  is  in  degrees. 


Since  Y  =  1.0  for  0C  =  10  deg,  the  empirical  constant  K  for  an  arbitrary  cone  angle 
0C  with  bluntness  ratio  0.167  (in  the  range  6  <  <  20,  5  <  0C  <  30  deg)  is  the 

product  of  Eqs.  (8)  and  (9). 
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Figure  15.  Correction  to  the  parameter  K  for  cone  angles  other 
than  10  deg. 


A  graphical  representation  for  several  finite  cone  angles  is  shown  in  Fig.  16.  It  should 
be  noted  that  this  relationship  is  still  limited  to  cone  nose  bluntness  values  of  0.167. 
As  in  all  empirical  approaches,  the  degree  of  confidence  decreases  as  one  departs  from 
the  original  physical  model.  Additional  data  are  needed  in  Fig.  16  before  great  confidence 
can  be  placed  on  the  results  of  the  analysis. 
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Figure  16.  Graphical  representation  of  the  parameter  K  at 
various  conejpigles  and  Mach  numbers. 

6.1.3  Estimate  of  the  Effect  of  Nose  Bluntness  on  the  Parameter  K 

By  using  the  present  experimental  data,  the  parameter  K  corrects  (in  part)  the  error 
introduced  by  assuming  a  model  controlled  by  sharp  cone  inviscid  flow-field  relationships. 
Since  these  data  were  obtained  on  a  blunt  cone  with  £  =  0.167,  the  previous  analysis 
could  not  expect  to  be  valid  for  cone  bluntness  at  values  other  than  0.167. 

Rather  tedious  hand  calculations  using  the  characteristics  solutions  shown  in  Fig.  10 
and  two  additional  solutions  for  a  25-deg  cone  suggest  that  the  variation  in  the  parameter 
(x)(Ps/P«,)  as  compared  to  a  £  =  0.167  solution,  behaves  as  A0"  where  n  is  on  the  order 
of  unity.  The  ratio  of  local  shock  angle  0S  for  a  blunt  cone  solution  to  that  of  a  sharp 
cone  value  was  calculated  and  normalized  by  the  value  at  £  =  0.167.  The  results  for  two 
cone  angles  (at  two  Mach  numbers)  are  shown  as  a  function  of  £  in  Fig.  17.  For  the 
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10-deg  cone  solution  the  results  show  a  strong  function  of  £  and  a  weak  function  of 
Mach  number  in  the  range  from  14  to  20.  The  solutions  for  the  25-deg  cone  were  only 
available  for  values  of  £  5>  0.23,  However,  at  this  value,  (0s)b/(0s)s  was  very  near  unity, 
and  it  was  assumed  that  the  value  at  £  =  0.167  would  also  be  unity.  The  numerical  effort 
and  time  to  develop  the  results  shown  in  Fig.  17  exceeded  the  effort  of  all  previous 
results  discussed  in  Section  6.  Because  of  the  problems  indicated  in  this  figure,  the  effect 
of  model  bluntness  discussed  later  was  restricted  to 

0C  -  10  deg 

0  <  £  <  0.3 

It  was  also  assumed  that  the  results  shown  in  Fig.  17  could  be  extended  from 
=  14  to  =  6.5,  and  the  parameter  K  was  adjusted  from  the  previous  value  by  multiplying 
by  the  relationship,  shown  in  Fig.  17, 


Z  =  f  (8 


(10) 


2  ■  Km J^Km!^, 


Sym  ^9  Mg) 


Figure  17.  Estimate  of  the  effect  of  nose  bluntness  on  the  parameter  K 
for  10-  and  25-deg  cones. 
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This  empirical  parameter  can  be  expressed, 
for  0  <  $  ^  0.13, 

Z  »  1.3  -  0.26291  £  -  5.3788  £2  -  48.418  £3  -  248.23 
for  0.13  ^  i  <  0.175, 


Z  =  1.0 


(ID 

(12) 


and  for  0.175  ^  $  <  0.30, 

Z  ~  0.52849  +  2.7257  £  (13) 

As  is  probably  evident  in  the  previous  discussion,  the  analysis  is  basically  a  rather  large 
extrapolation  of  a  meager  amount  of  data  to  conditions  beyond  the  combination  of  model 
geometry  and  flow  conditions  investigated  experimentally.  A  straightforward  analysis  using 
blunt  body  method  of  characteristics  solutions  combined  with  additional  experimental  data 
could  clarify  the  situation. 

6.2  SLIDING  FRICTION  DRAG 

The  mechanism  of  the  forces  induced  because  of  sliding  friction  is  beyond  the  scope 
of  the  present  analysis.  However,  it  was  suggested  that  the  analysis  at  least  include  this 
effect  qualitatively,  so  the  term  Cdf  in  Eq.  (1)  was  retained  and  defined  as 

CDp  =  FD/q^A  (14) 

where  FD  is  an  arbitrary  friction  drag  expressed  in  lbf  and  can  be  set  at  any  value. 

6.3  PRESSURE  DRAG 

Hopefully,  the  largest  component  of  drag  in  Eq.  (1)  would  be  inviscid  pressure  drag, 
Cdp-  In  hypersonic  flow  with  sphere  cones,  Newtonian  theory  has  proven  useful  in 
determining  this  parameter.  A  closed  form  solution  of  axial  force  is  presented  in  Ref. 
13,  and  for  a  =  0  it  can  be  expressed  as 

CDp  =  Cy-—  [(l  -  -y-  COS20c^  (2  sin20c)  +  £2  COS20cJ  (15) 

where  Cpmax  for  Mach  numbers  greater  than  6  has  the  value  (7  +  3)/(y  +  1).  The  ratio 
of  specific  heats  for  air  (1.4)  was  used  in  all  calculations  of  the  present  analysis.  There 
would  be  no  pressure  drag  acting  on  the  cylinder  portion  of  the  model  at  a  =  0. 
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6.4  SKIN  FRICTION  DRAG 


At  the  velocities,  range  pressures,  and  model  sizes  proposed  for  a  full-scale  track 
system  it  is  probably  correct  to  assume  that  the  boundary  layer  will  be  fully  turbulent 
over  a  major  portion  of  the  model.  For  a  sharp  cone,  the  turbulent  skin  friction  drag 
referenced  to  the  base  area  can  be  expressed  by  (Ref.  14) 


where 


'D  V 


Q.0776F  /VA1  8  /pe\°  8  /T-\°  58  /He\°-58 

(Refc)0-2  yTe  J  \H*j 


cot  Qc 


(16) 


-§f-  =  0.5  +  0.5  (Tw/T J  (TJTe)  +  0.0374  Me*  (17) 

The  skin  friction  coefficient  decreases  with  increasing  nose  radius.  However,  mass 
addition  caused  by  ablation  greatly  reduces  the  skin  friction  as  determined  through  the 
mass  addition  term,  F,  in  the  above  equation.  Typical  variations  of  F  with  enthalpy  for 
typical  materials  are  shown  in  Fig.  18.  When  Eq.  (16)  was  checked  against  more  recent 
turbulent  boundary-layer  solutions  available  in  VKF  for  several  cone  geometries  and  flow 
conditions,  values  of  CDy  from  Eq.  (16)  from  25  to  40  percent  above  the  VKF  solutions 


Figure  18.  Variation  of  mass  addition  term  with  enthalpy  for 
several  ablative  materials. 
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resulted.  Equation  (16)  does  not  estimate  the  skin  friction  component  of  drag  on  the 
cylinder  portion  of  the  model.  Since  this  apparent  "error"  in  Eq.  (16)  is  in  the  positive 
direction,  this  contribution  to  the  total  skin  friction  drag  would  at  least  be  approximated 
by  use  of  Eq.  (16)  without  modification.  It  should  also  be  noted  that  the  term  Cdv 
is  usually  quite  small  in  relation  to  Cdp  or  CDt,  so  a  25-  or  40-percent  error  in  CDy 
would  not  be  significant.  The  stagnation  enthalpy  term  H/RT0  plotted  in  Fig.  18  should 
be  evaluated  using  real  gas  relationships  to  determine  the  parameter  F  in  Eq.  (16).  Its 
effect  is  always  to  reduce  the  contribution  of  skin  friction  to  total  vehicle  drag. 

6.5  BASE  PRESSURE  DRAG 

For  drag  forces  defined  as  in  the  present  case,  the  contribution  of  a  pressure 
differential  in  the  model  base  region  to  the  free-stream  value  produces  a  drag  component 
given  by 


Cdb  "  W"  O  '  Pb/Pj  (18) 

There  are  numerous  data  correlations  for  the  ratio  of  Pb/p„,-  For  hypersonic  turbulent 
boundary-layer  conditions  in  the  base  region  of  sharp  or  slightly  blunted  slender  cones, 
the  correlation  shown  in  Fig.  19  can  be  used.  The  subscript  e  refers  to  local  inviscid 


Figure  19.  Base  pressure  data  correlation. 
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edge  conditions;  therefore,  the  ratio  Pb/p„  must  be  obtained  at  a  given  local  Mach  number 
by  multiplying  Pb/Pe  from  Fig.  19  by  a  value  of  Pc/'P^  obtained  from  an  inviscid  flow-field 
solution  for  a  sharp  cone  or  a  characteristics  solution  for  a  blunt  cone  at  the  end  of 
the  conical  portion.  The  present  method  uses  an  empirical  curve  fit  for  both  of  these 
methods  in  the  present  range  of  interest.  It  should  be  noted  that  Eq.  (18)  neglects  any 
effect  of  the  expansion  around  the  corner  of  the  cone  to  the  cylindrical  portion  of  the 
model  and  also  neglects  the  influence  of  the  track  on  the  average  pressure  in  the  base 
region.  This  problem  is  amenable  to  experimental  study.  For  certain  large  cone  angles, 
Eq.  (18)  will  predict  a  negative  base  drag  component  (thrust).  This  occurs,  however,  only 
when  forebody  pressure  drag  is  quite  large. 

6.6  INDUCED  PRESSURE  DRAG 

At  turbulent  boundary-layer  conditions,  the  effect  of  viscous-induccd  pressure  drag 
on  total  vehicle  drag  is  usually  less  than  one  percent.  For  this  reason,  induced  pressure 
drag  is  neglected  in  the  present  analysis. 

7.0  RESULTS  OF  THE  ANALYSIS  AND  CONCLUSIONS 

The  relationships  derived  and  discussed  in  Section  6.0  were  programmed,  with  the 
necessary  auxiliary  relationships,  for  the  AEDC-IBM  360/370  computer.  Parametric  studies 
of  the  different  variables  were  then  conducted.  Input  information  required,  and 
corresponding  empirical  limits,  were  as  follows.  All  calculations  were  performed  for  a  range 
ambient  temperature  of  300°K. 

Model  velocity,  fps,  7,410  <  <  28,475 

Range  pressure,  torr,  no  limits 

Track  diameter,  in.,  no  limits 

Track  rail  width,  in.,  no  limits 

Cone  half-angle,  deg,  5  <  0C  <  30 

i  0  <  £  <  0.3  for  0C  =  10  deg 

t  |  =  0.167  for  Qc  =£  10  deg 

Model  weight,  Ibm,  no  limits 

Friction  drag,  lbf,  no  limits 

Model  wall  temperature,  °K,  no  limits 

Mass  fraction,  0  <  F  <  1.0 
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With  this  information,  the  program  calculates  the  necessary  free-stream  conditions 
and  each  component  of  drag  in  Eq.  (1).  Also,  the  ballistic  coefficient  is  calculated  using 
the  expression 

B  =  W/CdtA  (19) 

An  extremely  large  number  of  combinations  of  input  variables  could  be  devised  from 
the  information  listed  above.  Only  a  few  were  examined,  and  the  results  are  discussed 
below. 

The  effect  of  model  wall  temperature  on  the  term  CD  v  [Eq.  (16)]  and  total  drag, 
CpT  [Eq.  (1)],  was  studied  for  a  10-deg  cone  at  model  velocity  of  20,000  fps.  The  effect 
of  varying  wall  temperature  from  300  to  3000°K  on  the  term  Cdsi/CDt  was  only  about 
2  percent.  All  subsequent  calculations  were  therefore  performed  for  a  wall  temperature 
of  811#K  (1000°F).  A  similar  small  effect  of  the  mass  fraction  parameter,  F[Eq.  (16)], 
was  also  noted,  and  subsequent  calculations  were  performed  for  a  nonablating  model  (F 
=  1.0). 

The  majority  of  the  remaining  calculations  to  be  discussed  were  performed  for  a 
10-deg,  slightly  blunted  cone  with  %  =  0.167.  This  model  corresponds  to  the  available 
experimental  data  and  results  in  the  greatest  confidence  in  the  parameter  K  (Section  6.0). 
Track  diameter,  TD,  was  arbitrarily  chosen  for  most  solutions  at  1 0.0  in.  and  rail  width, 
RW,  was  scaled  directly  from  the  present  experimental  apparatus  (Fig.  1)  with  a  resulting 
value  of  2.308  in.  Unless  otherwise  indicated,  all  solutions  assume  zero  friction  between 
model  and  rail. 

The  contribution  of  shock  interaction  drag,  expressed  as  CDsi/CDt,  as  a  function 
of  model  velocity  and  range  pressure  is  shown  in  Fig.  20.  A  very  strong  influence  of 
model  velocity  and  a  weaker  dependence  on  range  pressure  is  shown  for  this  model  and 
track  configuration.  Increase  in  total  drag  due  to  model  bow  shock  interaction  varies  from 
about  9  percent  to  as  much  as  30  percent  of  the  total  deceleration  force  of  the  model 
over  the  range  of  velocities  and  pressures  shown.  The  same  information  is  shown  in  a 
different  format  in  Fig.  21.  In  order  to  calculate  the  ballistic  coefficient,  B,  model  weight 
is  required.  All  solutions  shown  in  Fig.  21  use  a  model  weight  of  65  lbm.  This  represents 
the  present  estimated  launcher  limit  for  a  lO.O-in.-diam,  two-stage  light  gas  gun  at  a  launch 
velocity  of  20,000  fps.  Although  the  lower  velocities  shown  in  Fig.  21  could  have  used 
a  higher  model  mass,  the  solutions  retained  a  value  of  65  lbm  to  prevent  the  introduction 
of  another  variable.  Subsequent  calculations  for  launchers  of  diameters  other  than  1 0  in. 
and  launch  velocities  other  than  20,000  fps  often  use  the  maximum  model  weight  for 
a  given  combination  of  diameter  and  velocity  calculated  from  the  empirical  relationship 

Wmax  =  65  (4£)3  (20,000)2/(VJ2  (20) 
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where  Wmax  is  in  pounds  mass,  TD  is  in  inches,  and  is  in  feet  per  second.  The  ballistic 
coefficient  varied  from  as  low  as  775  to  as  high  as  970'  lbm/ft2  over  the  range  of  range 
pressures  and  model  velocities  shown.  The  nonlinear  influence  of  model  velocity  should 
be  noted.  The  optimum  velocity  was  10,000  fps,  with  lower  ballistic  coefficients  resulting 
for  both  higher  and  lower  model  velocities.  This  was  probably  due  to  the  fact  that,  for 
this  model,  the  term  Cds,  was  close  to  the  term  CDy  in  Eq.  (1).  For  a  given  velocity, 
the  relative  effect  of  range  pressure  is  shown  to  decrease  as  it  increases  toward  one 
atmosphere. 

The  effect  of  model  bluntness  ratio  and  velocity  on  shock  interaction  drag  on  a 
10-deg  cone  is  shown  in  Fig.  22.  The  calculations  were  performed  for  a  constant  range 
pressure  of  100  torr.  A  very  large  influence  of  both  bluntness  ratio  and  velocity  can  be 
seen  to  exist  with  the  contribution  of  model  bow  shock  increasing  to  greater  than  40 
percent  of  the  total  drag  as  the  model  becomes  sharp  and  velocity  increases  to  24,000 
fps.  A  more  realistic  upper  velocity  would  be  20,000  fps.  The  effect  of  range  pressure 
and  model  bluntness  is  shown  in  Fig.  23  for  this  velocity.  Comparing  Figs.  22  and  23 


Figure  22.  Effect  of  bluntness  ratio  and  velocity  on  10-deg  cone 
shock  interaction  drag. 
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indicates  the  relative  influence  of  model  velocity  and  range  pressure  at  a  given  bluntness 
ratio.  It  should  be  noted  that  the  term  CDs[/CDt  becomes  smaller  as  £  increases  because 
of  a  large  increase  in  the  pressure  drag  term,  CDp,  in  Eq.  (I)  and  not  because  of  a  decrease 
in  shock  interaction  drag,  Cds,. 


ID  -  10.0  in. 
Pgj.torr  RW  •  2.308 in. 


Figure  23.  Effect  of  bluntness  ratio  and  range  pressure  on 
10-deg  cone  shock  interaction  drag. 


The  effect  of  track  diameter  on  shock  interaction  drag  for  10-deg  models  of  constant 
nose  radius  is  shown  in  Fig.  24.  These  solutions  were  run  for  a  model  velocity  of  20,000 
fps  and  a  range  pressure  of  100  torr.  It  was  necessary  to  vary  nose  bluntness  ratio,  £, 
from  0.036  to  0.250  to  produce  these  solutions.  Rail  width  was  increased  or  decreased 
in  direct  proportion  to  track  diameter  with  the  value  of  2.308  in.  being  retained  for  the 
10.0-in.  track  diameter.  It  is  seen  that,  for  a  given  nose  radius,  the  effect  of  increasing 
track  diameter  is  small,  but  for  a  given  track  diameter,  decreasing  nose  radius  results  in 
a  greater  contribution  of  shock  interaction  drag  as  a  percent  of  total  model  drag.  This 
latter  conclusion  was  also  obtained  from  Figs.  22  and  23  and  again  results  from  the  increase 
in  Cdp  rather  than  from  a  decrease  in  Cdsi  as  nose  radius  increases. 

All  previous  solutions  have  utilized  a  model  with  a  cone  angle  of  10  deg.  The  effect 
of  varying  cone  angle  and  model  velocity  on  shock  interaction  drag  in  terms  of  CDsi/CDt 
is  shown  in  Fig.  25.  All  of  these  solutions  were  performed  for  a  bluntness  ratio,  £,  of 
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Figure  25.  Effect  of  cone  angle  and  velocity  on  shock  interaction  drag. 
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0.167  and  a  track  diameter  of  10.0  in.,  with  a  rail  width  of  2.308  in.  Range  pressure 
was  held  at  100  torr.  A  most  unusual  behavior  with  increasing  cone  angle  is  seen  for 
several  velocity  values.  A  combination  of  high  velocity  and  cone  angle  results  in  shock 
interaction  drag  being  greater  than  40  percent  of  the  total  vehicle  drag.  Conversely,  for 
small  cone  angles,  the  effect  of  velocity  is  quite  small,  with  the  shock  interaction  drag 
component  being  only  about  14  to  18  percent  of  the  total  drag.  The  Mach  number  values 
listed  with  the  velocity  values  in  Fig.  25  are  for  a  free-stream  (range)  temperature  of 
300°K  as  well  as  a  free-stream  pressure  of  100  torr.  As  a  matter  of  interest,  range 
measurements  obtained  by  the  McDonnell  Douglas  Corporation  on  a  45-deg  sharp 
cone-cylinder  model  (Ref.  1)  are  also  shown  in  Fig.  25.  Although  these  data  are  in 
qualitative  agreement  with  the  present  analysis,  this  may  be  fortuitous  since  the  shock 
mechanism  is  completely  different  from  that  presently  assumed  (Figs.  7  and  9),  and  the 
track  data  may  have  included  rail  friction  forces,  which  are  not  considered  in  Fig.  25. 

The  procedure  of  directly  scaling  rail  width  to  track  diameter  from  the  values  shown 
in  Fig.  1  to  larger  track  diameters  was  somewhat  arbitary.  Other  design  considerations 
may  require  that  rail  width  be  scaled  by  a  factor  other  than  unity.  From  a  viewpoint 
of  shock  interaction  drag,  it  would  be  advantageous  if  the  scaling  constant  were  less  than 
unity.  This  is  shown  graphically  in  Fig.  26.  The  influence  of  rail  width  in  terms  of  the 
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Figure  26.  Effect  of  rail  width  on  maximum  ballistic 
coefficient  and  shock  interaction  drag. 
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parameter  CDsi/Cqt  and  ballistic  coefficient  are  both  shown.  The  calculations  were 
performed  for  a  10.0-in.  track  with  a  10-deg  cone  of  bluntness  ratio  0.167  at  a  velocity 
of  20,000  fps  and  range  pressure  of  100  torr.  Present  two-stage  light-gas  gun  development 
suggests  that  future  systems  of  10-in.  diameter  launching  a  model  at  20,000  fps  would 
be  limited  to  a  model  weight  no  greater  than  65  Ibm.  This  value  would  yield  the  greatest 
ballistic  coefficient  and  is  the  value  used  in  the  calculation  shown  in  Fig.  26.  The  parameter 
4(RW)/7r(TD)  (ratio  of  total  rail  width  to  track  circumference)  is  also  indicated  in  Fig. 
26.  This  parameter  is  utilized  later  in  graphical  and  tabular  presentations.  A  very  large 
reduction  in  the  shock  interaction  drag  in  terms  of  total  drag  and  available  ballistic 
coefficient  is  observed  as  rail  width  is  decreased. 

All  previous  solutions  have  been  obtained  with  the  assumption  of  zero  sliding  friction 
between  the  cylinder  portion  of  the  model  and  the  track  rail.  This  was  done  because 
the  present  experimental  program  was  designed  only  to  study  shock  interaction  and  no 
accurate  estimate  has  as  yet  been  made  on  the  sliding  friction  term  (Cdf)  in  Eq.  (1). 
Track  curvature,  nonsymmetric  ablation,  model  weight,  and  materials  all  influence  friction 
as  well  as  model  velocity.  From  Eq.  (1),  total  "track-related"  drag  force  would  be  CdSj 
plus  Cdf-  This  parameter  in  terms  of  CpT  and  the  previous  parameter  Cdsi/Cdt  is  shown 
as  a  function  of  absolute  track  friction  drag  in  Fig.  27,  The  solutions  were  obtained  at 
velocities  of  8,000  and  20,000  fps,  with  other  input  information  tabulated  in  the  figure. 
Since  the  friction  force  is  included  in  CDt  [Eq.  ( 1 )] ,  the  parameter  CDsi/CDt  decreases 
with  increasing  track  friction  drag.  The  parameter  (CDs,  +  CDf)/CDt  is  seen  to  increase 
much  faster  for  the  lower  velocity  case  than  for  the  higher  velocity  value.  For  an  arbitrary 


Figure  27.  Effect  of  friction  drag  at  two  velocities. 
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track  friction  of  2000  lbf  at  a  velocity  of  8000  fps,  the  total  track-related  deceleration 
force  represents  about  76  percent  of  the  total  vehicle  drag  for  the  model  configuration 
indicated  in  Fig.  27.  Although  this  is  a  rather  severe  penalty,  the  actual  friction  force 
would  probably  be  much  less  than  2000  pounds  at  this  velocity  and  the  penalty  quickly 
becomes  smaller  as  friction  is  decreased.  A  more  realistic  estimate  of  the  relative  magnitude 
of  shock  interaction  and  track  friction  drag  is  shown  in  Fig.  28.  For  lower  velocities 
and  friction,  friction-related  drag  is  small  compared  to  shock  interaction  drag  and  increases 
to  about  the  same  magnitude  at  higher  velocities  and  absolute  friction  forces.  The  extreme 
slope  of  the  parameter  Cp  F/Cp  ^  in  Fig.  28  suggests  the  need  for  better  definition  of 
these  expected  track  friction  forces.  All  additional  calculations  discussed  herein  assume 
zero  friction. 


Figure  28.  Ratio  of  friction  to  shock  interaction  drag  as  a 
function  of  velocity. 
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Planning  plots  for  four  typical  model  geometries  and  range  conditions  are  shown 
in  Fig.  29.  The  solutions  are  shown  as  ballistic  coefficients  as  a  function  of  model  weight 
and  track  diameter.  Also  shown  is  the  present  estimated  limit  for  launch  weights  of 
two-stage  light-gas  guns.  With  this  as  the  maximum  available  ballistic  coefficient,  the 
summary  plot  shown  in  Fig.  30  was  prepared.  The  penalty  for  shock  interaction  drag 
ranged  from  about  14  to  43  percent  for  these  cases.  The  advantage  of  increasing  track 
diameter  in  obtaining  greater  ballistic  coefficients  can  be  seen.  It  should  also  be 
remembered  that  track  friction  forces  are  assumed  to  be  zero  for  these  calculations. 

A  slightly  different  type  of  planning  plot  is  shown  in  Fig.  31.  To  minimize  the  total 
number  of  variables,  all  these  solutions  were  obtained  for  a  model  velocity  of  18,000 
fps,  but  model  weights  were  calculated  for  launch  velocities  of  17,000  and  20,000  fps 
over  the  launcher  diameter  (track  diameter)  range  indicated  on  the  figure.  As  can  be  seen 
in  Fig.  25,  the  error  involved  in  using  a  model  weight  calculated  for  a  launch  velocity 
of  17,000  fps  and  then  "flown"  at  18,000  fps  (for  a  5-  or  10-deg  cone)  is  not  large. 
Model  cone  angles  of  5  and  10  deg  represent  typical  full-scale  reentry  vehicle  (RV)  values 
which  might  be  fired  in  a  full-scale  ablation  range  facility.  As  expected,  available  ballistic 
coefficient  is  a  strong  function  of  track  diameter  and  a  weaker  function  of  rail  width. 
Maximum  available  ballistic  coefficient  varies  from  about  500  to  greater  than  4000. 

■  LO  deg 
£  ■  0.167 

.  Pq,  ■  lOOtorr 


a.  9C  =  10  deg,  V.  =  20,000  fps 
Figure  29.  Planning  plot  for  four  typical  model  geometries 
and  range  conditions. 
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d.  9C  =  1U  deg,  V 
Figure  29.  C 
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The  two  extreme  values  (minimum  and  maximum  B)  in  Fig.  31  were  chosen  to 
perform  trajectory  calculations,*  and  the  results  are  shown  in  Fig.  32.  These  trajectory 
solutions  utilize  an  aerodynamic  drag  curve  as  a  function  of  velocity  calculated  from  Eq. 
(1).  and  in  addition,  Eq.  (1)  was  solved  with  the  term  CDsr  set  equal  to  zero  (no  shock 
interaction).  The  trajectory  solutions  were  run  for  a  time  interval  sufficient  to  allow 


Figure  31.  Planning  case  with  zero  friction  forces. 


a  full-scale  range  distance  of  50,000  ft  to  be  traversed.  The  resulting  ballistic  coefficient 
as  a  function  of  range  length  is  shown  in  Fig.  32a,  and  model  velocities  are  shown  in 
Fig.  32b.  For  these  models,  track  dimension,  and  range  conditions,  the  penalty  for  utilizing 
a  constraining  track  is  no  greater  than  13  percent  in  terms  of  ballistic  coefficient  and 
1300  fps  in  terms  of  reduced  model  velocity.  These  solutions  were  run  assuming  no  sliding 
friction  forces  were  present  between  the  model  cylinder  and  rail  surface. 


♦Trajectory  solutions  were  obtained  by  J.  P.  Billingsley,  VKF/ADP. 
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5-deg  Cone  -  •(Fig.  31) 
Launch  Velocity  *  17, 000  fps 
TD  •  16  in. 

Model  Weight  -  368.5  lb 
4IRW)hr(TD)  *  0.10 
i  *0.167 
Pm’  lOOtorr 


10-deq«Cone  -  t/  (Fig.  31) 
Launch  Velocity  •  20t  000  fps 
TD  •  6  in. 

Model  Weight  *14.0  lb 
4(RW)br(TD)  *0.30 
$  *0.167 
Pd, "  100  torr 


C^-f  (Vgj,  Model,  Track  Dimension) 

from  Eq.  (1)  for  All  Solutions 

_ i  -  -J _ I _ I _ I 

0  10  20  30  40  50 

Range  Length,  ft  x  10'; 


a.  Ballistic  coefficient 
Figure  32.  Results  of  trajectory  solutions. 


Finally,  an  additional  set  of  solutions  is  shown  in  tabular  format  in  Table  5.  Initial 
conditions  and  table  nomenclature  are  as  follows: 

Input  Data 

V„  =  UINF  =  12,000,  16,000,  20,000,  and  24,000  fps 
=  PINF  =  10,  100,  and  1000  torr 
0C  =  THETAC  =  5.0,  7.5,  and  10.0  deg 

%  =  RN/RB  =  0.167  for  0C  =  5.0  and  7.5  deg 

=  0,  0.10,  0.167,  0.2,  and  0.3  for  0C  =  10  deg 

4(RW)/(TD)  =  RAIL/CIR  =  0.1,  0.2,  0.3,  and  0.4 
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Cds,/CDt  =  CDSI/CDT 


b.  Veloi 
Figure  32. 


Table  5.  Shock  Interaction  Drag  for  Aeroballistic  Track  Systems 


TntlAC 


HAIL/CIR 


CUAE.RU 


CDS! /COT 


5.0 
s.o  _ 
5.0 
5.0 
5.0 


12U{)0.0  10.00 

12VU0.0  10. 00 

12000.0  10.00 

120O0.0  10.00 

160O0.0  10.00 

16j.it). 0 _ 10.00 

16000.0  10.00 

16UU0.Q  10.00 

200OU.0  10.00 

20000.0 _  .  10.00 

20000.0  10.00 

>OJJO.U  _ _ 10-UQ 

24000.0  10.00 

24000.0 _ Ik.  00 

2*000.0  lit. 00 

1200.0x0 _ 100x114— 

12000. 0  100.00 

iytiiio.0  lpu.au 

120O0.0  100. 00 

16000L.U  .  .  100-00 

16000.0  100.00 

.16000.0  _ _ Lilk.00_ 

16000.0  100.00 


20000.0  100.00 

20000. 0..  _100.xkU _ 

20000.0  100.00 

24000. 0  100x00 _ 

24000.0  100.00 


24000.0  100.00  5.0 

12000. 0 _ lOOkxOU _ 5x0 

12000.0  1000.00  5.0 

__120uO»U  .  1UOO.00 _ S.iL 

12000.0  100O.O0  5.0 


16000.0  1000.00  5.0 

16OO0.U  1000. 00  6.0 

16000.0  loOO.OO  5.0 

20ooo.o  lOoo.oo  5..0. 

200J0.0  1000.00  5.U 


200O0.O  1000.00 
24000.0  1000.UO 
240U0.0  1000.00 
24000.0  1000.00 
24000.0  1000.00 


0.167 


0.167 


0.167 

0.167 

0.167 


0.167 

_ 0x167 

0.16/ 

—D.IOjL 

0.167 


0.167 


0.167 

_ 0.167 

0.16/ 


0.167 


0.16-/ 

..  0.1  &Z _ 

0.16/ 

_.__0xL6/- _ 

0.16/ 


0.16/ 

_ 0x16.7 _ 

0.167 


0.1.0 

_  o.io _ 

0.20 


0.16/ 


o.oaor 

_ kxilfliU- 

0.0607 

u.nabT 

0.0744 

_ 11.0/44 

0.0744 

_ U.OIAA- 

0.0716 

D.0716 

0.0716 

_ 0.0716 

0.0702 

_ 0x07.02- 

0.0702 

_ 0.0698 

0.0698 

_ 0.U6VB 

0.0698 

0.0639 

0.06*59 

_ 0x0*39  _ 

0.0639 


0.0612 


0.0612 

0.0599 

0.0599 


0.0599 

_ M.062S- 

0 . 0629 
_  0.0629 

0.0629 
0.05/2 


0.0572 

_ 0x05.7  2_ 

0.05/2 

0.0547 

0.0547 


U.U547 

_ 01x0.534 . 

0.0534 
_  0.0534 

0.0534 


0.0038 


0.0114 


0.0845 


0.0921 
Q.09S9 
0.0035  0.0779 

0.0071 _ Qx.0B.lft 

0.0106  0.0850 

-0...0 1 4 1 _ 0..  0Bb5_ 

0.0039  0.0754 

0 .0.011 _  0.0/93 

0.0116  0.0831 

U. 1)154 _ Q.087U 

0.0044  0.0747 

0.01 34 _ 0.0836 

0.0178  0.0880 

_UxOU38 _ 0x0.7-36  . 

U.0076  0.0774 


0.0449 


0.0152 

0.0035 

0.0071 

0.0106 

0.0141 


U.0077 


0.0154 

_ 0 . 0044 

0.0089 


0.0178 

0-0038 

U.0U76 

_ 0.0114 

U • 0 1 52 
U.0035 


0. 007 1 
0.0106 
0.0141 

_  0x0038- 

0.0077 


0.0154 

O.OU44 

0.0089 

.0x0.1.34- 

0.0178 


0.0850 
0 . 0674 
0.0709 

0.0745 _ 

0.0780 


0.0689 


0.0766 

-0.0643 _ 

0.0688 


0.0777 
0.0667  _ 
0.0705 

0.0743 _ 

0.0781 

.0' 


0.0643 
0.06/8 
0.0714 
0.0585  _ 
0.0624 


0.u7ul 

0.0578 

0.0623 

0.066/ 

0.0712 


0.1236 

0.1583 

0.0453 

11.0867 

0.1247 

_U.1S9.6l. 

0.0511 
-JO. 0974  . 
0.1392 
0.1 77S  . 
0.0596 
0.1597- 
0.2022 

_ 0x0515  - 

0.0981 
0. 14112 
0.1787 
-0.0524 
0.0996 
0.1423 
0.1812 


0.1120 


0.2015 

_0x(169 1_ 

0.1294 


0.2292 

0.0569 

0.1077 

0.1532 

0.1944 


0.1099 
0 .1562 
0.1980 
0 .0659 
(K1238 


0.2203 

0.0770 

0.1430 

0.2002 

0.2502 


Table  5.  Continued 


UINF 

- EliS _ 

PlNF 

— hwh _ 

THE T AC 
— Ufcfi. _ 

RN/HB 

HAIL/C1H 

CDAERO 

CUS1 

COT 

CD5I/CUT 

12U00.0 

12000.0 

10.00 
_ 10.00 

7.5 

7.5 

0.167 

n.11,7 

0.10 

(1.20 

0.0904 

0.0984 

U.U042 

0.1026 

-0*1069 _ 

0.11 11 

0.0414 
.  0.0795 

0.1147 

0.1474 

0.0532 

_ O-lull _ 

12000.0 

_ 12000.0 

10.00 

- HUM _ 

7.5 
_ 7.S 

0#  167 
_ 0 *.1*1 

0.30 

9.40 

0.09d4 

0.0904 

0.0127 

0.0170 

16000.0 

16UU0.D 

lo.oo 
- - 10-00 

7.5 

7.5 

0.167 

0.167 

0.10 

_ 0.20 

0.0944 
_ 0.0944 

0.0053 
_ 0.0106 

0.0997 

o.ioso _ 

U) 

to 


16000.0 


20000.0 

20000.0 


20000.0 

20000.0 


24000.0 

_2400(U0_ 

24000.0 

P4UOQ.O 


12000.0 

12000.0 


12000.0 


-ill.  DO _ 


10.00 

lo.QQ 


-Zx5_ 


10.00 

10.00 


10.00 

HUM 


7.5 
.  7x5 
7.5 
7.5 


-0.167 


0.161 

0.167 


¥*40 


0.167 

0.147 


0.10 

0.20 


10.00 

-10.00 


100.00 

10Q.QU 


0.30 

0.40 


-0.09-44 _ 

0.0931 

0.0931 


u.0159 

JU0.212. 


0.10 

U.2Q 


0.0931 

U.0931 


100,00 


7.5 

7.5 


7.5 


0.167 

-0x167 


0.167 


0.30 

9.AJ)  _ 

0.10 

0.20 


0.0929 

0.0920 


0.0060 

O._0.13fi_ 

0.0204 

Q.U272 


0.1103 

0.1156 


0.0999 

-0.1067 


0.30 


0,0929 
0.0929 
0.0069 
0 . 0069 


0.0085 

0.0171 


0.1135 

0.1203 


0,0069 


0.0256 

U.0341 

0.0042 

0.0065 


0.0127 


0.1014 
0.1099 
0.1185 
JUL27.0  . 
040911 
-0.0954 


0 • 1444 
0.1837 
0.0600 
-0.1274- 
0.1796 
0.2260 


0.0996 


0.0841 

- 0_.  1553. 

0.2161 
-0.2609- 
0.0466 
_ 0-0091 


_ 12000.0 

_ 100.00 

_ Zjl5 

0.167 

0.40 

0 1 0809 

..  31  70 _ 

0*1 Q19 

16000.0 

100.00 

7.5 

0.167 

0.10 

0.0824 

0.0053 

0.0078 

0.0605 

16000.0 

_ 100. 00 _ 

_ 7.S 

—0.15/ 

0.20 

0.0874 

0.01U6 

0. 1)931 

n  *  1 1  a  i 

16000.0 

100,00 

7.5 

0.167 

0.30 

0,0024 

0.9159 

.0.0904 

0.1619 

_ 16000. 0 

1  Oii  .(10 

_ 7.5 

-0.167 _ 

0.4U 

0.0824 

0.0212 

0.1037 

0.2048 

20000.0 

100,00 

7.5 

0.167 

0.10 

O.08QB 

9. 0060 

0.0076 

0,0776 

_ 20000.0 

_ 100.00 

_ 1.5 

0.167 

0.20 

0.0800 

0.0136 

0*0944 

20000.0 

10U.00 

7.5 

0.167 

0.30 

0.0800 

0.92U4 

9.1012 

0,2015 

_ 20000.0 

_ 100.00 

.7*5 

0.167 

0.40 

0.0808 

0x0272 

0*1080 

D.Klf) 

24000. 0 

100,00 

7.5 

0.167 

0.10 

0.0802 

0.0085 

0.0807 

0.0962 

_ 24000.0 

- 100. 00 _ 

7.5 

0.167 

0.21) 

0.0002 

0.0171 

0*0972 

0.1756 

24000.0 

loo.oo 

7.5 

0.167 

0.30 

0.0002 

0.025b 

0.1050 

0.2421 

_ 24000.0 

_ IQU.OO _ 

7.5 

0.167 

0.40 

0.0802 

0.0341 

0. 1 143 

0.9Q87 

12000.0 

1000.00 

7.5 

0.167 

0.10 

0,0796 

0.00*2 

0.0039 

0,0507 

_ 12000.0 

1000.00 

_ lib 

0.167 

0.20 

0 ■ 0  796 

U.OOtfS 

0.0881 

0.0965 

12000.0 

1000.00 

7.5 

0.167 

0.30 

0.0796 

0.9127 

0.1)924 

U. 1300 

- 120QD.0 

1 OOU. UP 

_ 7.5 

0-167 

V-40 

0.0790 

U.0170 

0.0906 

0.1 760 

16000.0 

1000, 00 

7.5 

0.167 

0.10 

0,0749 

9.0953 

0.0002 

0.0661 

_ 16000.0 

-1000.00 

_ 7.5 

0.167 

.  0.20 

0.0749 

U.OlUb 

0.00^5 

16QU0.0 

1000. OU 

7.5 

0.167 

6.30 

0.0749 

9.9159 

0.9900 

0,1753 

_ 16000.0 

1000.00 

7.5 

0.167 

6.40 

9,0249 

0.0962 

0 .9900 

20000.0 

1000.00 

7.5 

0.167 

0.10 

0.0730 

0.9960 

0.9790 

0.0851 

_ 20000.0 

1000.00 

_ 7.5 

6.20 

0.0730 

U.  1)130 

0.0866 

o.l 4?n 

20000.0 

iooo.oo 

7.5 

0.167 

U.30 

0.0730 

0.0204 

0.0934 

0.2183 

_ 20.0  0  (U  0 

looo.oo 

_ 7.5 

0.167 

O.A0 

0 .0730 

0.0272 

0.1002 

0.7713 

24000.0 

1000.00 

7.5 

0.167 

0.10 

O.0721 

0.0005 

0.0007 

0,1058 

.  _ 24QU0.U 

—1009.00 

7x5 _ 

0.  *7 

6.20 

9.9221 

O.U171 

0.08*2 

24000.0 

1000.00 

7.5 

0.167 

0.30 

9.0721 

0.0256 

0.0977 

0.2620 

- 24000.0 

iooo.oo 

_ 1x5 _ 

0-167 

- iUM _ 

-  0.0721 _ 

- 0.6341 _ 

0.1063 

0.3713 
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Table  5.  Continued 


C6 

co 


UINF 

FPS- 


12000.0 

120U0.O 

12UUO.O 

12UUQ.U 

16000.0 


PiNF 

IQHtf 


1 ME  I  AC 
OLlim _ 


PN/HB 


HAIL/CIR 


CUAEHO 


COSI. 


CUT 


CDS 1 /COT 


10. 00 
.10.00 
10.00 
10.00 
10.00 
»o-»u 


10. 0 
ie..u_ 
10.0 
10.0- 
lo.o 
lu.o 


0.0 

-O.u- 


o.o 

_o.u_ 

o.o 

ft  -  0 


U.10 

U«|>0 

0.30 

o.*n 


O.09B9 

-0*0989- 

0.0969 

0.0989 


0.10 


0.096S 

U-U965 


0.0061 
n.n  imp 
0.0242 
-U.0323- 
0.0111 
11*0222 


0. 1070 
O.llSi- 


0.1232 

0.1313 

0.1076 

»-» IHH 


0.075S 

0*3-406- 

0.1968 

0.2*02 


0.1035 
0  - 1 »7* 


16U00.U 
16000 • 0 
2(1000. 0 
20000.0 
200U0.0 

*fl_ 


10.00 
10.00 
10. UO 

. lo.oo. 
lo.  uo 
10- no 


10.0 

lo.o 

10.0 

_ 10.0- 

10. 0 
lu.o 


6.30 

0.0965 

0.0334 

0.1299 

0.2572 

U  .40 

Q-OQhfi 

U>Q44fa _ 

0.1411 

_ 0.3159. 

0.10 

0.0962 

0.0148 

0.111U 

0.1332 

U.20 

(1.0962 

(1.0296 

0.1258 _ 

Q  .2352. 

0.30 

0.0962 

0.0444 

0.1406- 

0.3157 

Ua40  _ 

0.0962 

_ 0-0592 _ 

_ 0.1554 _ 

0.3B09 

24000.0 

10.00 

10.0- 

0.0 

0.10 

0.0967 

f>*OUO.U 

in^ao 

10*0 

0.0 

6.20 

_ 0*09.61. 

24000.0 

lo.oo 

10.0 

0.0 

0.30 

0.0967 

24l)U0f0 

iU  +  U(J 

lu.u 

0.0 

-  .  -  0.40 

O*09b.7_ 

12000.0 

lOU.OO 

10.0 

0.0 

0.  10 

0 .0864 

_ I2uu0*fl - 

ioo.no _ 

lo.O - 

.  0.0 - 

0-0864 

0.U169 

0.0370 


0.0568 

0.0757 


0.0061 

0.0162 


-0.114 
0.1535 
-0.1724 


0.0945 

0.1026 


0.1636 
.201.4 
0.3700 
JJL.-43.92- 


0.0655 

U..157S 


12000.0 
12000.0 
16000.0 
16000. 0 
16000.0 
1600U.U 


100.00 

lQU*ilO- 

100.00 

luu.oU 

lOO.CO 

100-00 


10.0 
_La. 
10.0 
lU 
10.0 
10-0 


0.0 

-0.0- 

0.0 

n.Q 


0.30 

u.*n 


0.0864 
0 . 086. 


0.0 

n  -  u 


0.10 

_a.2o. 

0.30 

U.«n 


0.0630 

0.0630 


0.0830 


0.0242 
-D..I 
0.0111 
(U0222 
0.0334 
0-0446 


0.1107 
0.1188 
0.0942 
L.053- 
0.1164 
0.1276 


0.2190 

0.2722 

0.1182 

Q.2116 

0.2870 

0-4*93 


20000. 0 
20OU0.0  . 
20000.0 
20000*0 
240UO.U 

2*0oo.o 


100.00 

_ 10-0. 00 

100.00 

_ 100,00.. 

100.00 

luu-ou 


10. 0 


0.0 


0.10 

6.20 


O.JO 

-4L.4U- 


0.10 

0.20 


0.0820 

0.0148 

0.0968 

0.1528 

0.08>0 

U.U296 

0.1116 

0.2652 

0.0820 

0.0444 

0,1264 

0.3511 

0.0820 

0.0592 

fl.1412 

0.4192 _ 

0.0619 

0.0189 

0.1008 

0.1877 

0 . OHIO _ 

0.0379 _ 

0*1197 _ 

_ 6.3162 _ 

24000.0 
2400Q„.0  - 
12000.0 
.12000.0. 
12000. 0 


100. uo 

loo.^o.o 

1000,00 

...1000,00 

100U.OO 


10.0 

-10.0 

10.0 

10.0 


0.0 

-fl-.il 


0.30 

0..Q 


0.0 

_0.A_ 

0.0 


0.10 

0.20 


0.0819 

_a*0B19_ 

0.0785 

0.0785 


0.0568 

0.0757 


0.30 


0.0765 


O.Q0B1 

JU0.l62_ 

0.024? 


0.1386 

0.1S76 


0.0866 

0.0947 


0.1028 


0.4095 

JU.48Q5- 

0.0932 

0.1707 


0.2358 


16000.0 

16OO0.U 

16000.0 

1BU00.0 

20000.0 

loou.oo 
1000, 00  \  . 
1000,00 
ItiOO.OU 
looo.ou 

ll)«0 

|dfO 

— 

0.0 

0.0 

6.10 

0.20 

0.0745 

0.0111 

0.0223 

0.0657 

0.0968 

0.1300 

0.2301 

10.0 
10.0  . 
10.0 

0*0 

u.O 

0.30 

U.40 

0.0745 

0 .0745 

0.0334 

0.0446 

0*1080 
0.1191 _ 

0.3095 

2*3742 _ 

0.0 

0.10 

0.0730 

0.0148 

0.0878 

0.1683 

_ 20000.0 _ 

1060. 00 _ 

10.0 

0.0 

0.2u _ 

0.0730 

U-0296 _ 

0.1026 

0.2883 

2ouo0.o 

200o0.u 

24UU0.0 

2400U.ll 

24000.0 

24OJU.0 


1000.00 

1000.00 

1000.00 

1000,00 

luOd.00 

lfltfV.UO 


10.0 

10.0 

10.0 

10.0 

10.0 

lOitf 


0.0 

o.n 


0.0 
_fl  «_U_ 


0.  Jo 
0_*40_ 
0.10 
6.20 


0.0730 
710_ 
0.0725 
0.0725 


0.U592 


0.6189 

0.0379 


0.0914 

0.1104 


0.0 

0.0 


0.30 

V«40 


0.0725 
(i . 0725 


0.0568 

A.H757 


0.1293 
Oil  Mi. 


0.2069 

0.3*31 


4392 

5109 
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12000.0 

_12000.0 _ 

12000.0 
12000.0 
16000.0 
IftOOO-O 
16000.0 
16000.0 
20000.0 
20000.0 
20000.0 
20000-0 
24000.0 
24000.0 
24000.0 
24000.0 
12000.0 
12000.0 
12000.0 
12000.0 
16000.0 
16000.0 
16000.0 
MAitA  n 


20000.0 
20000.0 
20000. 0 

-26000.0 _ 

240U0.0 


10.00 


10.00 


10. 00 


10.00 
—10.00- 
10,00 
ly.oo 
100.00 
.  100.00 
100.00 
100.00 
100.00 
.  loo .  0JL_ 
100.00 
i  ml  a  n 


100. 00 
100.00 
100.00 
lOu.oo 
lOo.oo 


TmETaC 

llfcfi. 

10.0 

10.0 _ 

10.0 


10.0 

llt.O _ 

10.0 

_L0_.  0 _ 

10.0 


KAlL/ClK 


0.100 

-0.100 

0.100 

_(U100_ 

0.100 

0-100 

0.100 

o.lou 

0.100 

0.100 

0.100 

0.100 

0.100 

-0...1PJUL 

0.100 

0.100. 

0.100 

o.lou 

0.100 

0*_1Q.<L_ 

0.100 

0.100 

0.100 


0.100 

0.100 

0.100 


160O0.0 
16000.0 
16000.0 
20000.0 
20000.0 
'lAAllrt  n 


20000.0 


24000. 0 
j>40ll0.0 
24000.0 


1000.00 

1000.00 

1000.00 


1000.00 


1000,00 


1000. oo 
1000.00 
1000.00 


0.100 

o.lou 

0.100 


0.30 

Jl.4(L- 

0.10 


0.4U 

.  -JL.3Q _ 

0.10 


cuaeko 

oTToTi” 

0.1075 
0,1075 
0.1075 
0.1050 
0.1050 
0.1050 
-  0. 1030 
0.1047 
-- UjIQA-Z- 
0,1047 

0.1047 
0.10S2 
0 .1052 
0.1052 
-0-.1Q52 
0.0050 
0.0950 
0.0950 
0.0950 
0.0916 
0.0916 
0.0916 


0.0905 
0.0905 
0.0905 
_0»09ji4  . 
0.0904 


TT’A'l'l.l 


U.0904 


10. 0 

_1J)..  0 _ 

10.0 


III 


0.0871 

0.0871 


0.0871 


0.0831 

0.0B31 

0.0831 

0.0616 

0.0816 


0.0816 


0.0810 

0.UB10 

0.0810 


0.0072 

_ 0.0144 

0.0216 
_ 0.0288- 
0.0101 

_ 0.02O2 

0.0J03 

_ 0.04-tL4_ 

0.0136 

_ u.0271 

0.0407 

_ 0.0543 

0.0174 

_ 0.0-3.49- 

0.0698 

0.0524 

0,0072 

_ 0-.0144 

0,0216 
.  0.0288 
0.0101 
_P.0-2.OiL- 
0.0303 


0.0271 


0.0349 


U.069B 
0.0072 


0.0144 
0,0216 


0,0288 


U.U202 

_ t.0303_ 

0.0404 

0.0136 

0.0271 


0.0543 

_ 0.0.17.4 

0.0349 

_ 0.0524 

0.0698 


0.1147 
Pt 1219 
0.1291 
0.1363 
0.1151 
0.1253 
0.1354 

-  Otl.455 _ 

0.1183 
0.1319 
0.1454 
0.1590 
0.1227 
0.1401 
0.1751 
0.1576 
0.1022 
.-Oil  09.4 
0.1166 
0.1238 
0.1017 
0.1118 
0.1219 


0.1177 

0.1448 

0-t_l  0  7  8_ 
0.1253 


0.1602 

0.0943 


0.101S 

0.1087 


0.1159 


0.1033 

0.1134 

0.1235 

0.0951 

0.1087 


0.1358 

0.0985 

0.1160 

0.1334 

0.15U9 


COSI/CIIT 

0.0627 

_ 0.11 80 

0.1671 

0.211.1 _ 

0.0877 

_ 0-1614 

0.2240 

_ 0,2-780 . 

0 a  1146 

0.2057 _ 

0.2798 

_ 0.3413 

0.1422 

-0.2491 

0.3989 

_ P..3323  __ 

0.0704 

_ 0.1315 

0.1851 

.  -0..2325 _ 

0.0994 

-9,.  18_09 _ 

0.2468 


0.2306 

: _ 0.3100 

0.3.747 

0.1618 _ 

0.2786 


0.4358 

0.0 


0.1418 
0.1965 


0.2483 


0.1957 

_ 0,2674  _ 

0.3274 

0.1425 

0.2495 


0.3994 

—0.17-72 

0.3011 

0.3925 _ 

0.4629 
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Ui  HF 

PPQ 

HI‘*F 

ll)MD 

ThElAC 

nl-fci 

RN/H8 

HAIL/CIR 

COAt.HU 

CObI 

CUT 

CDS I /CUT 

12000.0 

10.00 

10.0 

0.167 

0.10 

0.1231 

0.0061 

0.1291 

0.0470 

1  7  0  Dli.i) 

1  U  .  Ill) 

10.0 

0.167 

0.70 _ 

0 • 1231—  . 

..  0.0121  . 

0.1352 

0.0698 

12000. 0 

10.00 

10.0 

0.167 

0.30 

0,1231 

0.0162 

0. 141 J 

0.1289 

] ^UUQ.Q 

lu  t  hu 

1Q.Q 

n.  167 

o.«o 

0.1231  _ 

.  0.0243 

.0.147J 

0.1648 

16000. 0 

10.00 

10.0 

0.167 

0.10 

0.1205 

U.0085 

0.1291 

0 . 0659 

_ 16000. Q _ 

_ 10-00 _ 

10.0 

"-167 

_ 6*26 _ 

_ 0-17U5 _ 

_ U.U17U _ 

».137t. _ 

_ 11.1237 

16000.0 

10.00 

10.0 

0.167 

O.JO 

0.1205 

0.0255 

0.1461 

0.1747 

1 6000. 0 

lu.00  _ 16*6 

0.167 

U.«0 

0*1205.  . 

.  0.0340 

0.1546 

0.2202 

20000.0 

10.00 

10.0 

0.167 

0.10 

U  •  1202 

0.0114 

0.1316 

0.0865 

200b0.ll 

1U.00 

io. u _ _ 

0.167 

0.20 _ 

0.1207 _ 

_ 0*6226.  ... 

0*1430 

0.1594 

20U00.0 

10.00 

10.0 

0.167 

0.30 

0.1202 

0.0342 

0.1544 

0.2214 

_ ?umlil.Q _ 

_ 10-00 

_ 10.0 _ 

_ 0*16/ _ 

_ 6*46 _ 

_ 6*1262 _ 

0-0456 _ 

0*  lbHfft 

0.775u 

iyi 

<-6 


24000.0 

74000.0 


240U0.0 

2401)0.0 


12000.0 

1701)11-0 


12000.0 

12000. 0 


16000.0 

160U0.P 


16000.0 

160U0.Q 


10.00 

qp.uu- 


10.0 

10.0- 


10.00 


100.00 

IQO-OU 


10.0 
-12* 
10.0 
1U.U 


167 

167 

167 

16/_ 


10 
20 
JO 
4  fl¬ 


it? 

162- 


100,00 
100.00 
100.00 
100. Ou 


10.0 


167 


10 

2IL. 


12o7 

12llZ_ 

12o7 

1207. 

1105 

11 115 


0.0147 
.0293- 
0.0440 
.9*9586 
0.0061 
0.0121 


0 

0 

0 

.9* 

0 

-0. 


1394 
1500  . 
1647 
2134.. 
1166 
1222. 


0.1062 

0.1954 

0.2670 

0.3269 

U.QS20 

n.noon 


30 


1105 


0.0162 


1266 
0.134B 


0.1414 


100.00 

10,0,.  00 


10.0 

0.167 

0.10 

0.1071 

u.0085 

0.1156 

0.0736 

10. U 

6*167 

_ 0.20 

0.1071 

0.0LZ6  ... 

..  0.1241  .. 

...  11. 13.7 1— 

10.0 

0.16/ 

0.30 

0.1071 

0.0255 

0.1326 

0.1924 

10.0 _ 

_ 0-167 _ 

_ 6*46 _ 

0.1071 

11.0340 

_ 6*1411 _ 

_ 0.2511 

20000. 0 

70000. 0 


20000.0 

70000.  II 


24000. 0 
2*000-0 


100.00 

1 00. UU 


10.0 

10.0 


167 

1 6/_ 


10U.00 

loo. ou 


10.0 

lo.o 


167 

162. 


luO.OO 

ioo.no 


10.0 

m.o 


167 

162. 


10 

20_ 

30 

4fl_ 


1060 

1660. 

1060 

1060- 


0.0114 

0.0226 


1174 

1266 


0.0 J42 
0.0451) 


10 

20 


1059 

1059 


0.0147 

0-0203 


1402 

1516 

1205 

1362 


0.0970 

0.1770 

0.2436 

-0.3007— 


0.1215 

0.2169 


24000.0 
-24000. 
12000. 0 
12000.0 
12000.0 
1701)0. D 


loo.oo 

lQU.OU 

loou. oo 

JJ 
loou.oo 
1000-00 


10.0 

m.OL 

10.0 


10.0 

m.o 


167 

IU167-. 


30 

46. 


1059 

1059 


0.0440 

u.0566 


167 

162. 


10 

20- 


1026 

1626- 


11.0061 

11.0121 


1499 

.1645 

1067 

1146, 


167 

162. 


30 

46- 


1026 

1626 


u.0162 

u.0243 


1209 

1269- 


0.2934 
.3564- 
0.0556 
0.1656- 
0.1507 
0.1913 


16000.0 

1000. 00 

10.0 

0.16/ 

0.10 

0.0986 

0.0065 

0.1071 

0.0794 

lbUUO.O 

1000*00 

10.0 

0.167 

0.70 

u  ■  6966 

0.01Z6 _ 

_ 0*1156 

0.1472 

16000.0 

1000.00 

10.0 

0.16? 

0.30 

0.0986 

0.0255 

0.1241 

0.2056 

16000. 0 

1000.00 

1  Da  u 

(1*16/ 

0.40 

0.0966 

0.0340 

0.1326 

0*2566 _ 

20000.0 

1000.00 

10.0 

0.16/ 

0.10 

0.09/1 

0.0114 

0.1065 

0.1050 

_ 20000.0 _ 

1000.00 _ 

_ 16*6 _ 

_ 6*16/ _ 

_ 6*26 _ 

_ P.U97J _ 

0.0225 _ 

0.1199 

0-1902 _ 

20000.0 

20000.0 


24000.0 

24000.0 

24000.0 

7401)0.0 


lOUU.OO 

lOuO.OO 


lOOU.UU 
1000. 00 
lOUU.OO 
1660. 00 


10.0 

-16t6 _ 

10.0 
ASL, 
10.0 
10.0 


O.lo/ 

0.30 

0.0971 

0.0342 

0.1313 

0.2604 

U.40 

U.0971 

0.0456 

0.1427 

0  >3196 

0.167 

0.10 

0.0965 

U.0U7 

0.1112 

0.1318 

0.167 

0.20 

0.0965 

0.0293 

0.1259 

-0  *233.0. 

0.167 

0.30 

0*0905 

0.0440 

0.1405 

0.3130 

31*16/ - 

_ 0.40 _ 

_ 9*11965 _ 

_ 9*9656 _ 

_ 9*1662 _ 

0.3779 - 
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On 


uli.F 

_ FKS 

►  INF 

- Ulna _ 

IritlAC 
_ llfc.ll. _ 

hn/kb 

KA1L/CIR 

COAERO 

cost 

CUT 

COS 1 /COT 

■  120UO.O 

lo.oo 

10.0 

0.200 

0.10 

0.1337 

0.0060 

0.1397 

0.0429 

12000.0 

lo.uo 

--1U.U- 

... _ _ 

_ 0.20 

0.1 1^7 

a.oipa 

(1.1447 

0.0R24 

12000.U 

10.00 

10.0 

0.200 

0.30 

0.1337 

0.0180 

0.1517 

0.1186 

12000. 0 

10.00 

10.0 

__  0*20U _ 

_ U*AO _ 

_ 0.1337 

0.0240 _ 

0.1577 

0.1422 

16000. 0 

lo.oo 

10. 0 

0.200 

0.10 

0.1312 

0.0084 

0.1396 

0.0598 

- lliilflfl.il 

_ 10.00 

_ 10-u 

_ 0.200 

U.20 

(1.1.412 

A.OlliT 

0.1474 

16000.0 

lo.oo 

10.0 

0.2UO 

0.30 

0.1312 

0.0251 

0.1563 

0.1804 

loOUO.O 

10.00 

10.0 

-fl.2Q.u__ 

_  .0.40 _ 

_ u.1312 

_ IU1L33& 

0.1646 

0.2030 

2OUUO.0 

lo.oo 

10. 0 

0.200 

u.  10 

0.1309 

0.0111 

0.1420 

0.0785 

20000.0 

.  .1.0.00 

_lU.O_. 

_ 0  *21)0. 

_ 0.20 

0*1309 

0.0223 

0.1542 

0. 1457 

20000.0 

lu.00 

10.0 

0.200 

0.30 

0.1309 

0.0335 

0.1643 

0.2036 

0.9^41 

_ 2UUUO.O 

_ lu-uo _ 

_ ia«Q 

_ 0.20U 

U.4fl 

0 . 1309 

0.044A 

0.1755 

240UO.O  " 

lo.oo 

10.0 

0.200 

0.10 

0.1313 

0.0143 

0.1457 

0.0982 

24000.0 

_  -10.no 

_ ln.40 _ 

_  0.200. 

U.2U 

_ 0.1313 

0.0286 

0.1600 

0.1790 

24OU0.U 

10.00 

10. 0 

0.200 

0.30 

0.1313 

0.0429 

0.1743 

0.2464 

24000.0 

_  In.  oo  .. 

10.0 

_ 0*20 0_. 

0.1.40 

0.1313 

0.0573 

0.1804 

0.303? 

12000. 0 

100.00 

10.0 

0.200 

U.lO 

0.1212 

0.0060 

0.1272 

0.0472 

_ 0.0901 _ 

1200(1.0 

■liiu.oo.. 

10.0 

o.Pau 

_ 0.20 

_ 0.1212 _ 

_ u-oiyp 

_ 0*1342 

i2uuo.o 

160U0.U 

16000.0 

16000.0 

ibunn.u 


lOu.OU 


100, on 


lou.oo 

mumi 


—10.  <L 
10. 0 

_ IOjO.  . 

lo.o 

lo»o 


0.1212 

-0.1212- 

0.'U77 

0.1177 


0.0160 

0.0240 


0.1392 

Q.14S2- 


0.1177 


0.00B4 

0411107. 

0.0334 


0.1261 

0.1345 

0.1512 


0.1293 
0.1653^ 
0.0662 
1243 
0.2211 


20000.0 

100.00 

lo.O 

0.200 

0.10 

0.1167 

0.0111 

0.1278 

0.0872 

20000. 0_. 

_  100.00  . 

.  _in...o_ 

_ _ Q.2U0 

0*20  . 

0.1167 

0.0223 

0.1390 

0.1605 

20000.0 

100.00 

10.0 

0.2UU 

0.30 

0.1167 

0.0335 

0.1501 

0.2229 

20000*0  . 

—  1.017.00  . 

-lo.  0_ 

0*2017 _ 

-  0*45 _ 

0.1167 

0.0446 

0.1613 

0.2767 

24VU0.0 

100.00 

13.0 

0.200 

0.10 

0.1165 

0.0143 

0.1308 

0.1094 

_ 24(10  0.0 

_ IQU.tlO 

_ 10.0 

_ 0*2UU 
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Table  5.  Concluded 
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These  solutions  will  allow  the  interested  reader  to  make  graphical  studies  not  included 
in  the  present  report. 

The  present  experimental  and  analytic  investigations  indicate  the  approximate 
magnitude  of  aerodynamic  drag  which  will  be  induced  because  of  the  bow  shock  interaction 
between  the  track  and  a  model  traveling  down  a  constraining  track  system.  Although  it 
is  clear  that  additional  experimental  data  are  required  to  confirm  or  improve  the  analysis, 
the  following  parameters  (in  decreasing  order)  have  been  shown  to  be  important': 

1.  Model  Cone  Angle  and  Bluntness' 

The  relative  importance  of  shock  interaction  drag  in  terms  of  total 
drag  increases  with  increasing  cone  angle  and  decreasing  bluntness  ratio 
(Figs.  22,  23,  and  25)  except  for  small  and  6C. 

2.  Model  Velocity 

The  term  CdS|/Cdt  increases  very  rapidly  with  increasing  velocity 
(Fig.  20).  In  terms  of  practical  ballistic  coefficient,  an  optimum  velocity 
results  for  a  given  model  and  track  configuration  (Fig.  21).  A  combination 
of  high  cone  angle  and  model  velocity  results  in  greater  than  40  percent 
of  the  total  decleration  force  being  caused  by  model  bow  shock  interaction 
(Fig.  25). 

3.  Rail  Width 

The  contribution  of  shock  interference  drag  is  reduced  as  rail  width 
decreases  (Fig.  26). 

4.  Range  Pressure 

Shock  interference  drag  is  sensitive  to  range  pressures  at  levels  below 
about  100  torr  but  is  not  a  strong  influence  above  this  level  (Figs.  20  and 
21). 

5.  Track  Sliding  Friction 

Although  this  term  was  studied  only  qualitatively,  it  was  shown  to 
have  a  potentially  large  influence  if  the  absolute  value  is  greater  than  about 
100  Ibf  (Figs.  27  and  28). 
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NOMENCLATURE 

A  Reference  area,  7rrb  2 

B  Ballistic  coefficient  [ Eq .  ( 1 9) ] 

Ca;  Inviscid  axial-force  coefficient 

Ca  t  Total  measured  axial-force  coefficient 

Cd  b  Aerodynamic  base  drag  coefficient  [Eqs.  (1)  and  (18)1 

Cd  b  l  Baseline  drag  coefficient  measured  with  no  track  system  in  place 

Cp  f  Drag  coefficient  attributable  to  sliding  friction  between  model  cylinder  and  rail 
[Eqs.  (l)and  (14)] 

Cp  |p  Aerodynamic  viscous-induced  pressure  drag  coefficient  [Eq.  (1 )] 

Cd  p  Aerodynamic  pressure  drag  coefficient  [Eqs.  ( 1 )  and  (15)] 

CDsi  Drag  coefficient  attributable  to  shock  interaction  forces  between  model  and  rail 
[Eqs.  (1)  and  (3)] 

Cdt  Total  drag  coefficient  of  a  ballistic  range  model  traveling  in- a  constraining  track 
system  [Eq.  (1)] 

Cd  v  Aerodynamic  viscous  drag  coefficient  [  Eqs.  ( 1 )  and  (16)] 

CPmax  Maximum  pressure  coefficient  ■ 

Dsi  Shock  interaction  drag  [Eq.  (2)] 

F  Mass  fraction  term  in  Eq.  (16) 

FD  Sliding  friction  drag 

H*  Reference  enthalpy  defined  by  Eq.  (17) 

He  Gas  enthalpy  at  edge  of  boundary  layer 

Ho  Stagnation  enthalpy 

K  An  empirical  correction  factor  defined  by  Eqs.  (2)  and  (7) 
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L  Position  of  rail  leading  edge  in  relation  to  cone-cylinder  junction  (Fig.  1 ) 

Mi  ,£2  Characteristic  lengths  of  model  (Fig.  5) 

Cjj  Sharp  cone  length 

Me  Mach  number  at  edge  of  boundary  layer 

Mj  Local  Mach  number  after  first  reflection  (Fig.  9) 

M_  Free-stream  Mach  number 

Pb  Base  pressure 

pe  Local  pressure  at  edge  of  boundary  layer 

Po  Stagnation  chamber  pressure 

Po  Pitot  pressure 

Ps  Local  pressure  after  first  shock  reflection  (Fig.  9) 

pw  Wall  pressure 

p^  Free-stream  static  pressure 

q0  Stagnation  point  Fay-Riddell  heat-transfer  rate 

qL  Free-stream  dynamic  pressure 

R  Gas  constant 

Re  Unit  free-stream  Reynolds  number  per  foot 

Reg  Free-stream  Reynolds  number  based  on  characteristic  length  8 

Refc  Free-stream  Reynolds  number  based  on  sharp  cone  length, 

RW  Rail  width 

rb  Model  base  radius,  assumed  to  be  equal  to  track  inside  radius 

r„  Model  nose  radius 

Te  Temperature  at  edge  of  boundary  layer 
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T0  Stagnation  temperature 

Tw  Wall  temperature 

Free-stream  static  temperature 
TD  Track  diameter 

Ve  Velocity  at  edge  of  boundary  layer 

Free-stream  velocity 
W  Model  weight 

x  Cone  distance  influenced  by  shock  reflection  (Fig.  9) 

xi  Axial  distance  measured  from  nose 

Xb  Cone  distance  influenced  by  blunt  body  shock  reflection  (Fig.  10) 

xs  Cone  distance  influenced  by  sharp  body  shock  reflection  (Fig.  1 0) 

Y  Correction  to  the  parameter  K  for  cone  angles  other  than  10  deg[Eq.  (9)] 

y  Distance  from  cone  centerline 

Z  Correction  to  the  parameter  K  for  bluntness  ratios  other  than  0.167  [Eqs.  (10) 

through  (13)] 

a  Angle  of  attack 

y  Ratio  of  specific  heats 

6C  Cone  half-angle  (Fig.  9) 

0S  Local  model  bow  shock  angle  (Fig.  9) 

9's  Model  bow  shock  angle  after  first  reflection  from  rail  (Fig.  9) 

£  Nose  bluntness  ratio,  rn/rb 

p _  Free-stream  density 
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